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ABSTRACT 






A new method for design of control systems in the time domain is 
presented. The adjustable parameters, the roots and the zeros of the 
system are evaluated simultaneously either by hand or by computer 0 
The method is applicable to any system configuration and is able to 
handle a large number of adjustable parameters. The method is analy- 
tical, any degree of accuracy can be obtained. 
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PART I ANALYSIS 



CHAPTER I 

THEORETICAL CONSIDERATION OF ALGEBRAIC EQUATIONS 
1-1 . Introduction. 

The response of feedback control systems is governed by the complete 
set of the closed loop pole-zero configurations. Since the zeros are 
usually known, the only problem is the poles. For linear control systems, 
the problem of poles is conveyed to the roots of the characteristic equa- 
tion which is an algebraic equation. Several methods have been developed 

to this end such as Guilemin®s method*, root-locus method and Mitrovic’s 
2 3 4 

method * * . Guillemin°s method starts the synthesis by choosing a set 
of pole-zero configurations from the consideration of specifications and 
then finds out the network needed to give this chosen pole-zero configura- 
tion. The root-locus method is essentially a graphic method for one vari- 
able which has linear relation with the coefficients of the characteristic 
equation. Mitrovic’s method is also a graphic method but handles two 
variable coefficients. The method proposed in this paper is rather an 
analytic method and is able to handle any number of variables which are 
contained in the coefficients of the characteristic equation. The analy- 
sis is based on the root-coefficient relation of an algebraic equation, but 
a new technique of partitioning the equation is employed. 

1 - 2 . Root - coe f f ic ient Re lat ionship . 

For a n t * 1 order algebraic equation, there are n relations among the 
roots and the coefficients of the equation. If all the coefficients are 
actually known numerically, then all the roots are fixed in such a way that 



1 



' 





they satisfy the root-ccef f icient relations. Conversely* if a set of 
roots of the characteristic equation are chosen* then all the coefficients 
of the equation are fixed by those relations. However* if variables are 
contained in the coefficients* then the roots of the equation are varied 
in a definite relation to the variation of the variables. For example* 
if one coefficient of the characteristic equation is variable* then the 
relation between the roots and the variable coefficient can be plotted by 
root-locus method. But if this relationship is looked at from a different 
point of view* an analytical relation can be obtained. Consider a third 
order characteristic equation with one variable coefficient as shown in 
the following equation: 



s.s'-t uds + =o a-i) 

where £3^ and are known numerically* is a function of ~f) 

which can be varied. Assume -S, > -j> 2 anc * ~ are the three roots of 
the equation* then the root-coefficient relations are as follows: 



5, + -f % = Bi 



( 1 - 2 - 1 ) 



S,s, + + S,5, = f, (fil (1-2-2) 



S, S, S, = B„ 



(1-2-3) 



Since is adjustable* then equation (1-2-2) can always be satisfied 

for any choice of the roots by varying accordingly. However* 

equations (1-2-1) and (1-2-3) set constraints to the value of the roots. 
Namely* the three roots must satisfy the two equations (1-2-1) and (1-2-3) 
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If the three roots are treated as the variables of a set of simultaneous 
equations, then there are three unknowns but only two equations, consequent- 
ly one root is arbitrary. This arbitrary root can be chosen independently 
if there is no restriction to the value of f’i('fi) * a ^d therefore this 
root can be treated as an independent variable in the root-coefficient 
relations. In the root-locus method, the independent variable is /",(/?>) 
while in this point of view, the independent variable is conveyed to one 
of the roots of the characteristic equation. 

The above reasoning is based on the solution of simultaneous alge- 
braic equations, i.e., 5 / f 5 Z and are treated as independent 

variables. Since those variables are the roots of an algebraic equation, 
a complex root must be treated as two variables because of its conjugate 
root. If one variable of a pair of complex conjugated roots is expected 
to be treated as the independent variable, this variable cannot be chosen 
to be the root itself, but can be treated as follows. Assume - S f and -5^ 
are a pair of complex conjugate roots of the equation (1-1), those two 
roots can be transformed to other two variables by the following trans- 
formations : 

5 d-3) 

5,5^ = COn (1-4) 

Where ^ and iO n are known as the damping ratio and the natural frequency 
of the pair of complex conjugate roots —5, and . By this trans- 

formation, both ^ and tO n can be treated as independent variables. 

In equation (1-1), if f,{^) was also a fixed number, then all 
roots were fixed numbers by equations (1-2) and no root was arbitrary. 

The arbitrary roots arise from the introduction of variables to the coeffi- 
cients. If one variable is introduced to the coefficients, it has been 

3 



shown one root is arbitrary. For the same reason* if two variables are 

t h 

introduced to the coefficient of a n order characteristic equation* 
then two roots are arbitrary* and for YY] variables to the coefficients, 
]fY] roots are arbitrary. The root-locus method is the case of one vari- 
able while Mitrovic°s method is the case of two variables* but both methods 
are quite different ways of approach. It is because of this way of look- 
ing at the root-coef ficient relations* an analytic method can be derived 
and is capable of handling any number of variables contained in the co- 
efficients, 

1-3, Partition of a Characteristic Equation, 

If adjustable parameters are introduced to the coefficients of a 
characteristic equation* then this equation can be partitioned into two 
equations. One of them is formed by the arbitrary roots* the other is 
formed by the constrained roots. The constrained roots here are defined 
as those roots other than arbitrary. 

Consider a 5th order characteristic equation 

5 r -f- + B j S 3 +8 2 S i +6 ( S+ B 0 =o d-5-D 

Assume —5, — S? — S, - 5, and are ihe roots. The negative 

sign is for convenience. By expanding the equation 

(S + 5 ( )CS+5 i )Cit5 J J(i+3^)(5+5^J =0 (1-5-2) 

Set the corresponding coefficients of (1-5-1) and (1-5-2) equal s then the 
root-coefficient relations are obtained as follows; 

e 4 = S, + .S 1 + 5 J °' 6 ' l) 

B 3 = (5 itiKSjtiJ + SA-tV, +CV-S. + -S, °' 6 " 2) 
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* 






( 



& z =- J^(S,fS 2 ) + ( Sj+SJJ'Sj. i [(S,+Sj(5 J + ^)i SjS+JSj. (1-6-3) 

B, = >A 5 j^4 + [i S. + SjSjS^ + S^CSjtS^)] S s (1-6-4) 

Bo — 5,5^5^ (1-6-5) 

Assume the coefficients {3, and are adjustable, then equation (1-5-1) 

can be written as follows: 

S r + 84 . 5 % B35 3 4 4KA)-S* + fi (£,&)s -h 3 o ^0 (1-7) 

Where f (-#, 4) and ;f , ( # |x 4) are functions of f) k and which 

are adjustable parameters. f3 0 > Bj and £3^ are assumed given fixed 
values. A terminology, "controlled characteristic equation" is given to 
equation (1-7) on the reason that its roots are controllable by the ad- 
justable parameters. Since there are two adjustable parameters intro- 
duced to the coefficients, two roots can be chosen arbitrarily. Let those 
two arbitrarily chosen roots be ( -5>/ ) and ( ) , then the other three 

roots (-S,), ( -^ ) and ( -5^- ) are not arbitrary but constrained to the 
fixed coefficients ( g , and ) and the arbitrarily chosen roots 

as shown in the last section- By separating the arbitrary roots and the 
constrained roots, equations (1-6-1), (1-6-2) and (1-6-5) can be arranged 
as follows: 

From (1-6-1) 

3 j 4 3 4 + S 6 - = B 4 -(S, + SJ 0-8) 

From equation (1-6-2) 

B 3 = ( 3, +5 4 + S s )(S. + ^) + + (1_9) 
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( 



From equation (1-6-5) 



8 0 = (^sjc-V* 5 ^ (i-i°) 

It can be recognized that ( 5^ + \ + S s ) * ( Sj 5^ + ^ S^- 4 S j ) and 

( <> ) are the coefficients of a third order algebraic equation 

which has roots of -5j , and • Define 

C.= ^s*i <1 - 11 - 1) 

c,= S 3 S^ + S^S s -tS s S J (1-11-2) 

= + d-ii-3) 

« 

Then the equation of the constrained roots is formed by 

S 3 + C^S 2 + C, S + Co — O (M2) 



It must be noticed that equations (1-12) has the roots -S^ , -5^ and - 5^ 
Substitute the definitions of (1-11) into (1-8), (1-9) and (1-10) and 
manipulating, one obtains 



c. 



B U+s*) 



(1-12-1) 



C, =• Bj -Ci( s,+i)~ s,5 f 



( 1 - 12 - 2 ) 



C 



& /5,S* 



(1-12-3) 



Consider equations (1-12), > C/ and C 2 . are functions of the arbitrary 

roots and the fixed coefficients B„. 63 and 8 ^ of the controlled 

characteristic equation (1-7). Therefore the roots of equation (1-12) 
are also functions of the arbitrary roots - 5, and - 5 ^ and the fixed co- 
efficients. Assume there are no restrictions to the variable coefficients 

6 
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^ (f) 9 /,>nd f 2 (f) y j 7 ,)then t ^ ie controlled characteristic equation (1-7) 

is partitioned into the following equation 



(54S,)(s + s i )(5 , + c 2 s'+c,s +c 0 ) =o 

because both equations have the same roots. As Q q , C t a nd C 2 are 
functions of 5 , anc * S 2 , then the five roots of the original character- 
istic equation (1-7) are completely determined by equation (1-12). Equa- 
tion (1-12) is a third order algebraic equation and the coefficients are 
expressed analytically by (1-12-1), (1-12-2) and (1-12-3). Then the 
original 5th order characteristic equation is now reduced to a third order 
equation and by this reason, a terminology "REDUCED CHARACTERISTIC EQUATION 1 ; 
is defined for equation (1-12). 

In the above arguments it is assumed that the variables ^ and ^ 
are adjusted such that two arbitrary roots are chosen. Therefore, the 
variables ft), and f) ^ , which are the actual physical adjustable parameters 
of the system, must be determined. From equation (1-6-3), and (1-6-4) 
replace by f z ( ^ , ) and B, by f t ( ft, J $ J since they are 

assumed adjustable, and manipulating, one obtains 



r 

7r 




= C 0 + C,(S,+Si) +CjS,S z 



(1-14-1) 




C„( 5, + 3^) + C, S.S., 



(1-14-2) 



Substitute C s C, a nd C> from equations (l-l£) into (1-14), obtain 

0 I 



fzlMj - B^S.+S^W ^ S - +(S l + 5U i -25 l S e (V^) <1_15 " 1) 



f.tf.Aj = BUS.+SO/s,** 4 ^3 S ' 5 - - +5^) * 5,^U,Vj Z -(S,S 2 f(L-15-2) 

It can be seen from equations (1-15) that 'fj ( anc * fj ( ^ } are 

expressed analytically as a function of the arbitrary roots and the fixed 
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coefficients $ Q , B*. and 8^ • And ^ Iso can be expressed ex- 

plicitly as functions of the arbitrary roots. Therefore if a set of 
roots has been chosen from the consideration of the reduced characteristic 
equation, then the physical adjustable parameters are readily determined 
by equations (1-15) or (1-14). 

By this partitioning process it can be seen easily that the original 
fice root-coefficient relations of (1-6) are transformed into other five 
relations which are equations (1-12) and (1-14) or (1-13) and (1-15). 

For convenience (1-12) and (1-14) are repeated here: 

= i^-cs. + s*) 

C, — Bj-CiCs.+SJ-s.S* 

C 0 =■ B. As* o-w 

f, = C„ + C,(S +QS.S, 

f 2 — + C,S,S* 

In actual numerical computations, it is simpler to use equation (1-16) 
instead of (1-12) and (1-15), because of the progressive nature of the 
root-coefficient relations of algebraic equations. Namely, compute C 0 
and C t first, then using these calculated data to compute C f and then 

f, and • 

In general, the n root-coefficient relations of a nth order character- 
istic equation which has variables in coefficients can be arranged in any 
form as one wishes. The process derived here has certain significances 
in the design of control systems. First of all, the order of the character- 
istic equations which determines all the roots of the system is reduced by 
the number of the variable parameters. Secondary, the physical adjustable 
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parameters are expressed directly by the arbitrarily chosen roots* The 



stability and response of the system is determined completely by the re- 
duced characteristic equation while the adjustable parameters are deter- 
mined directly by the arbitrary roots. Since equations (1-16) are alge- 
braic equations and if no singularity is involved, then solutions of ft) f 



ways exist. Although this process does not guarantee the adjustable para- 
meter obtained in this procedure to be always realizable, the realizable 
region can be found. In the design of a control system, this is a logical 
procedure, namely the pole-zero configuration of the closed loop system is 
chosen first from the specifications, and then compute compensator para- 
meters directly from the chosen pole-zero configurations. For every 
choice of pole-zero configuration, a corresponding set of values of the 
compensator are obtained. This partitioning process effectively trans- 
forms the desired pole-zero configuration to be the independent variables 
of the design procedure. Moreover, most of the designs are based on a 
pair of dominant complex conjugate roots, then two of the arbitrary roots 
can be chosen (not necessary) as the dominant roots. The adjustable para- 
meters of the compensator therefore can be expressed explicitly as func- 
tions of the dominant roots if the dominant roots are defined. 

In this process, there are no limitations. For a control system of 
any order, any number of adjustable parameters and any type of compensa- 
tor; the same technique can be applied. 

1-4. General Procedures of Partitioning a Controlled Characteristic 



The example in the last section assumed that the adjustable parameters 
are contained in two coefficients. Now consider a cascade single section 



and which actually determine the physical adjustable parameters al- 




Equat ion. 
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R-C compensator as shown in the block diagram of Fig* 1-1, in which the 
compensator parameters enter all the coefficients* The characteristic 
equation of the compensated system is 

5 4 + ( B 2 + P)S } + ( B,+ 6 >p)S i +( + ^ = 0 (1-17) 



where B e = K 


, B, = <xb , B* = 


CLi b 










K 




V 




(5-hpJ 




5(5+a)(s+b) 

















Fig. 1-1 Single section cascade compensator* 

Ot = ? / , £ ; K » a and b are know numerically* 

Assume , - S 3 and - are the roots of (1-17). By expanding 

the equation 

(S+5 1 J(5 + S 2 )(S-*-5 3 )(5+4 ) .) = O (1-18) 

and setting the corresponding coefficients of (1-17) and (1-18) equal, 
one obtains the following four relations: 

B*+ P ~ -5, + +5 3 + -5+ (l-19-l) 



B,+ B a f> = (V^)(5 3 + S*) + S,S, + %% (1-19-2) 

B,p-+ = 5 A ( 5_> +-5^) -f 5 i ^.( S, + -%. ) (1-19-3) 

PBo — (1-19-4) 

Since there are assumed two adjustable parameters 0( and p of the compensa- 
tor, there are two arbitrary roots. Assume -5, 9 and - 5 2 are the 
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arbitrary roots, then ( + 5^. ) and ( 3^5^ ) are the coefficients 

of the reduced characteristic equation* Define 



C, = 3,4 S, 



( 1 - 20 - 1 ) 



C G = S 3 s + 



then the reduced characteristic equation is 

5 Z + C,S -hCo — O 

Substitute the definitions of (1-20) into (1-19), one obtains 



(1-20-2) 



( 1 - 21 ) 



1 


0 


- 1 


O 


(5,+s,) 


1 


-b* 


o 


5.5* 




-B, 


-Bo 


0 


3 , 5 , 




o 

J 



C, 

Co 

IP 

c< 



B, - 5,5, 
O 
O 



( 1 - 22 ) 



J 



Equations (1-22) are the partitioned four root-coefficients relations in 
which and C q define the reduced characteristic equation while p and X 
are the compensator parameters. These four dependent variables (C^, C q , 
p and ) can be solved in terms of the arbitrary roots and the fixed 

plant parameters ( g , g ( and (3^). For numerical computation purposes, 
only one dependent variable need be expressed explicitly as a function of 
the arbitrary roots, the others are computed readily by the previous cal- 
culations. Here C q is solved from (1-22) as shown in the following equa- 



tion 



where 



C„ = ir[ SA-(s,+s/- B, + 6 i (s,+s i )] 
A. — [ _ | 



(1-23-1) 



(1-23-2) 
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i 



















Now consider another example which introduces three adjustable 
parameters. Fig. 1-2 is a fourth order system with first and second 
derivative feedback. The parameters hj P h ^ and K (=s K^K 2 ) are assumed 




Fig. 1-2 Pure derivative feedback K s h^ 
and h^ are variables where 
K = K x K 2 

adjustable. The characteristic equation is 



^ + & 3 S s + ( + K z h 2 )s*+ ( 6 , 4 - K z h,) 5 + I < ~ O (1-24) 

where Bj = O. -f b + C , B>. , z= a_b + be + CO , and B, — CX.bC 

Assume the roots are ~3, s ~S Z » -S 3 ~ then the root coefficients 

relations are as follows: 



Bi= S, + + S 3 + 5* 


(1-25 -l) 


3* + K 2 h 2 r=r ( S, + 5 2 )S 3 + V* + ( S, + 5 3 ) 5* 


(1-25-2) 


B, + K 2 h 1 5, S z Sj -j- £ 5, + S 3 ( s, + 5 ^ } J 


(1-25-3) 


K = ^ s 3 S+ 


(1-25-4) 
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< 














Since there are three adjustable parameters (h^ h 2 and K) , three roots 
are arbitrary. Let these arbitrary roots be - S ( 9 » and - Sj . 

The reduced characteristic equation is of order one and its roots is - . 



Define S, = C , then die reduced characteristic equation is 
4 o 

5 K 0 =0 (1-26) 

Substitute the definition of S^ = into (1-25) , the root-coefficients 
relations become 

c„ = Bj- ( S,+ 5, + SjJ (1-27-1) 

K 2 h* = ( s,+ S,)Sj + S,5 i + C„(5 ,+5.+5 J )-B a (1-27-2) 

KA = i,S.S J + C.[S,S i + (S,+ 5 J )S 3 ] 0-27-3) 

K = 5,5,5 iG (1-27-4) 



Equations (1-27) are the partitioned root-coef f icient relations in which 
(1-27-1) determine the system response while the other three equations 
determine the adjustable parameters so that to give the roots determined 
by (1-27-1), 

From the three examples illustrated above, it can be seen that this 
partitioning process can be applied to systems of any order, any numbei of 
adjustable parameters and any system configurations. The important concept 
of this process is to separate the roots of the controlled characteristic 
equation into two parts: The arbitrary roots and the constrained roots. 

In the n root-coefficient relations, the arbitrary roots are treated as 
independent variables while the constrained roots and the adjustable para- 
meters are the dependent variables. The constrained roots are obtained 
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from the solution of the reduced characteristic equation while the ad- 
justable parameters are calculated from the arbitrary roots directly. 

The general procedures are summarized as follows: 

(1) Formulate the controlled characteristic equation*. Set 
the coefficients equal to the corresponding root-coefficient. 

(2) Define the arbitrary roots* and the coefficients of the 

reduced characteristic equation (C „ Substitute C^*.. into 

the equations in (1). 

(3) Solve for the adjustable parameters and C°s in terms of the 
arbitrary roots. 

1-5. Summary. 

This partitioning process effectively transforms the independent vari- 
ables from the adjustable parameters to the arbitrary roots of the character- 
istic equation. As the adjustable parameters (compensator) are determined 
from the desired closed loop pole-zero configuration, a design procedure 
for control systems based on the time domain can be derived. The parti- 
tioning procedures presented in the last section are applicable to all 
system configurations. For specific cases of a nth order system, general for- 
mulas can be derived. For the cases that some of the coefficients themselves 

*In discontinuous operation (2), the controlled characteristic equations 
may be different for different modes. As far as this step is concerned, 
each mode must be treated independently. 

+ The number of arbitrary roots is equal to the number of the independent 
adjustable parameters. The it dependent adjustable parameters may be less 
than the compensator parameters because of the specifications and other 
restrictions. For details, it is discussed in Chapter IV. 
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are treated as variables, general expressions from the partitioning 
process have been derived as shown in Table 1-1 to Table 1-6* For 
other cases, they are shown in the subsequent Chapters 0 In order to 
illustrate the applicability, the following examples are presented,, 
Example: Given a characteristic equation 

5^4 / 4? 4 1285*4 150OS + 3000 — O (1-28) 



which is unstable. Assume a network can be devised such that the co- 
efficients of the first, second and third derivative terms can be ad- 
justed, Then the controlled characteristic equation is 

4 \(o S* 4 fjS 5 ! f 2 S Z + 3000 =■ O (i_29 

where -f- 9 , f z and are variables. From Table 1-5, the coefficients 



of the reduced characteristic equation are: 

c, = e 4 -( S, + S;+5,J (1-30-1) 

C„= S„/s, 5, Si (1-30-2) 

The functions of the variable coefficients are: 

ft = C 0 (S ( 5 2 + ^5 3 45 3 Sj-f C,S,S 2 S s (1-30-3) 

+ + C.CS'St+^Sj+JjS,) + S.^Ss (1-30-4) 

f 3 = Co + C, (5,45 2 1-Sj) 4 ( S,5, +S;$ 3 i S 3 5,) (1-30-5) 



where (~S, ), ( ” 3 Z ) and ( - 5_j ) are the arbitrarily chosen roots. 
Since equation (1-29) is of the 5th order, one root must be real. Assume 
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this real root is ( -5^ ). For convenience let ( - S, ) and ( -5, ) be 
the dominant complex roots and express them by the following notations: 



5>, + S z — (1-31-1) 

.a- 

3,5i = n (1-31-2) 

Substitute (1-31-1), (1-31-2), = 16, and B> a = 3000 into the equa- 

tions (1-30), one obtains: 

C, — !(o — ( Cj n + S 3 ) (1-32-1) 

Cq— 300Oy/u)o (1-32-2) 

f,= C 0 (Ldn+2f(*)'}S i )+C l U)n S 3 (1-32-3) 

fz = C„(2fuJ n + S 3 ) + C,(u}*+ lfU) n S s ) 4 (1-32-4) 

f 3 — C 0 +C,(2j , u) rt fS }> ) + ( U)* i 2fUJ rt SjJ (1-32-5) 



The reduced characteristic equation is 

3 + Cl 5 "+ Co ::::::::: C> (1-33) 

4 

which determines all the roots of equation (1-29) and consequently deter- 
mines the stability boundary. By putting all the three arbitrary roots 
on the left -half plane, then from (1-3? 1) C q is always positive, the 
condition of stability for the arbitrary roots therefore from (1-33) is 

c, 

or f to ^ 2 ? 6 U n (1-33) 
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The inequality (1-33) implies that the three arbitrary chosen roots are 
confined in a certain region if the system is to be stable. This region 
is plotted in Fig. 1-3. Notice is the real part of the complex 

roots. If Sj and S 2 are complex the boundary is a line ( S = -8). If 
Sj and are real, the boundary is a point, (S = -8 +j O ). 

Assume the dominant root3 ( ~S, ) and ( -5, ) are desired such 
that all other roots must be to the left of the dominant roots, then 
the following inequalities must be satisfied 

fU> n (1-34-1) 

and C, > 2ftO n (1-34-2) 

Substitute the boundary conditions of (1-34-1) and (1-34-2) into (1-32-1) 
one obtains 

fu) n £: f (?/ 5 (1-35) 

Equation (1-35) defines the dominant root region which is also shown in 
Fig. 1-3. Within this region, choose = 0.5, C0 n = 4, = 4, then 

from (1-32-1) and (1-32-2) 

c, = a 

Co — 4-7 

The roots of (1-33) are: -4 + j 5.6 ( CO'^ = 6.85, J = 0.59) 

The root configuration for this choice of the arbitrary roots is also 
shown in Fig. 1-3. 

The next step of the design is to calculate the variable coefficients. 
Substitute ^ = 0.5, U) n =4, = 4, = 8 and C q = 47 into equations 

(1-32-3), (1-32-4) and (1-32-5), one obtains 

j- t = 2015 
- 696 
fj = 143 
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For another choice of the arbitrary roots* the coefficients of the 
reduced characteristic equation and the variable coefficients are com- 
puted in the same way* 

The above example assumed that three coefficients are variable. Other 
cases, such as a cascade compensator, feedback compensator, etc., can be 
carried out in the same procedure, since the expressions of the coeffici- 
ents of the reduced characteristic equation and the adjustable parameter 
are obtained readily. The stable root region air dominant oot region are 
important in design, the general methods for evaluating those two regions 
are discussed in the next chapter. 
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Fig. 1-3. Example for 5*+ loS^f fjS** f 2 -S^+ 13000 = 0 
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order 
of the 
equation 


coeff. of the reduced 
chara. equation 


functions of the 
variable coeff,. 


'fih 


c 2 = ( 

C, = Bj-C S,fSJ 

r _ Bo* 


fj = c 0 (;,+SJ +c,s,s z 
C e + C ( (5,iS 2 ) i C 2 S,S 2 


5 th 


<0=1 
Ci = 

f* Bo 

C °“ S,S* 


bth 


C, = 1 

C 3 = + 

C; = %-Cj (S,tSj-S,S* 

C ( = Bj, - C z (5, +5,) ~Cj S,S t 

r — JSk. 

“ 5>S* 


nth 


C r >-2 = « 

<-n-3 ^ i 

(3n-z “Cn-j (5,+^j-C,.^ 

Co-i* — S />-3 - C«-4 
• 
t 

c, = 6j- CAW) -c 5 ; f s 2 

r do 

° “ Vt 



Table 1-1 Functions of the variable coefficien 
functions contained in this table are assumed only B 

c Snd <“ S 2^ are 'the arbitrarily chosen roots and 
of the reduced characteristic equation. 

-f C,s i Co =c 



5 n + 






s f, and fi . The • 
and B 2 are adjustable. 
C's are the coefficients 
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order 
of the 
equation 


coeff. of the reduced 
cnara. equation 


Functions of the...... 

variable coeff. 


4 th 


^ = 1 ' , 
C, = O,- Cot 

C 0 — Bo / 5/^2 




% 

jnh 


c j = 1 

c*=v (s, + sj 

c, = [ B,- c.(s,«j]as. 
C 0 = B./SA 


Co+ C, (5,t 


(p th 


Of. = 1 

C } = B r - (5. + 5J 
Ci =■ - C* ( s,f .O ~ 

C, = [ &,-CoWSt)]/s,S* 


j3 = C,4 C,(5,t5j tCjJ^ 




Cr)~2 — I 

Cn-j — ~ C n-J- ( S,f5») 

C tj-4 — &r>-2 £ rt-3 C ” C|j.^ 5|5i 

* # » • 

t • * . 

\ • * 


/ 

t 




C 2 = B.-Cj(5,+iJ-C + S,S, 
c l = [ B,-C.(5,+S J )]/5,it- 
C„= B.A.Sj 


• 



Table 1-2. Functions of the variable coefficients f. and f* 

T ® f “ n ^ ions contained ^ this table are assumed only B 0 and B are” 

ar^thp ^ <\ S 2> are the arbitrarily chosen loots a^d C's 

re the coefficients of the reduced characteristic equation. 
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order 
of the 
equation 


coeff. of the redeced 
.cnara. equation 


Functions of tne 
variable coeff. 


4-t h 


C { - c B,-C o-S,S z ) As,*s 2 ) 

C 0 = Bo / 5, 5 Z 

• 






C,= l 

C.= B.-(S, + SJ 

C,= ( 6,-C„- C^Sj/CS.+Sz) 

C„= B./5.S. 


. f,= C. (*,«») + C,5,s. 
f 3 - C, + C 2 (S,+5,)4 Cj5,S 1 

1 

) 


(o +h 


qccl 

C-J - 6 s --(S,t5 z ) 

B*-C 3 (s l +S*)-S l S* 

C. = ( B t -C e - CtSA) fa St ) 
C? = 60 / 5,54 


nth' 


C„. a = 1 

C- n-s — Cft-i ( S,f5fc) 

=: ' Bp-x — O r >-3 ( 5, +^t) “ Cq-i ^,5* 

. • * 

* 

* 4 * * 

C 2 = B 4 -c j(V5 J-c 4 55 i 
C, = [ 

C» = B./s,S, 



Table 1-3 Functions of the variable coefficients ~ft and , 
The functions contained in this table are assumed only B 1 and B are 
adjustable. (-S^) and (-S~) are the arbitrarily chosen roots and C's 
are coefficients of the reduced characteristic equation. 



\ 
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order 
of the 
equation 


coeff. of the reduced 
chara. equation 


Functions of the 
variable coeff. 


3rd 


C, = i 

Co — B*-(S ,+ SO 


• 


4-th 


Co= B i - C, (5,t5^) - C 2 5,Si 
C, - 

Ci — 1 




5th 


Co = e t - ($,+*,) c, — 5,$ t c t 
C, = Sj - U+JJCi-UCj 
Cl — B*-(S,+ SJ 
Cj — / 


fo= c.s,s t 

ft = c.<s,-*s.)+c.s,s 2 




C-o — B*. — *-• #zC* 

C. = Bj — (Jtii) G — V>Cj 

• , • 




n+ h 


» • 

* « « • 

t 

« 






• 

C/).j ~ 

^-n-z — ) 


/’ , . 



and fj 



Table 1-4. Functions of the variable coefficients 
The functions in this table are assumed B and B, are adjustable 
simultaneously. Where (-S ) and (-S,) are the arbitrarily chosen 



v ° 1 -' ** uu are cne arbitrarily chosen 

roots and C's are the coefficients of the reduced characteristic equation. 
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order 
of the 
equation 


coeff. of the reduced 
chara. equation 


5 ih 


C,= 1 

c,= e*-u+V5,j 
c 0 - &*/$&$, 


6 th 


C } = 1 

Cl = 6 s-KS'+Sl+Sj) 

c, = a*-Cjt(s l fM5i)-c J (s,s I t«.vv.) 

C 0 

C 6 = Bp/s,5 t 5 3 


nth 


C ft _j — | 

e„..= b^. ,-C„(s,tMy 

n-ir — - Crt-4 (S|ti z i-Sj) ~C n _^ j^i 

n-(» = B„. 3 Cfj-s- (i, 

1 ' 

* , • • 


\ 


• • • • 

C, = B^-Qrs^tSj) -c^.s^y^,) 
C.= B./4,5 t S J 



functions of the 
variable coeff. 






o 

o 

-f 

r> 



Km 






A* 

v» 






vO 

$ 

or 






11 , ii 



n 

o 






<*> 



+ 

r> 



_v> 

VI 



*4" 

v/> 



_v» 

-WC- 



(I 

n 

0 

\A 

+ 

0" 

4 - 

Vy> 

u 

-h 

p 

v\ 

vY' 



Table 1-5. Functions of the variable coefficients f, f and -f 
The functions contained in this table are assumed B, , and B Z are 3 
adjustable simultaneously. Where (-S ), (-S ) and f-sj are t£e arbitrarily 

chosen roots and the C s are the coefficients of the reduced characterise 
equation. 



:ic 



\ 



i 
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order 
of the 
equation 


coeff. of the reduced 
chara. equation 


4+h 


+5 t 45 j) 

c, =. | 

• 


5- * 


S = 6j- c,U^U)-( s^+^+SjS,) 
S — B4 ( 5, + ^ + -5,^) 

Q = l 




c. = - G «,+ V*;) - Cj ( U iSJj tSjS.) 

c * =• %■- c, - c 4 

• « 

* H . • 


n 


t « • # 

^-fh4 = @0-1 ~ ( 

CfHj =| 

V 




- 



functions of tne 
variable coeff. 



II 

n 

9 



fSi 



& 



+ 

r> 

V/-> 

4 - 

& 



& 

V> 



vO 

„v> 

0 * 



S' 

o 

II II 



r> 

o 

_VA 

N^ 

+ 

•+• 



4 - 

o 

y\ 

t.V> 



p 

Vn 

£ 



T . a , Ta - b Jj 1_6 " Functions of the variable coefficients f 

!" n M i0n J C ° ntained ln this table are assumed B , B, , and B 
adjustable simultaneously. Notations are the same §s T lble 1-5 2 



f , and f z . 



and B2 are 
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CHAPTER II 



REDUCED CHARACTERISTIC EQUATION 

2-1. Introduction: The reduced characteristic equation determines all the 

roots of the controlled characteristic equation and consequently deter* 
mines the response of the system. Usually the zeros of the closed loop 
system are known, then the total performance of the system depends upon 
the choice of the poles (roots of the characteristic equation). In 
order to establish a guide for the choice of the roots* the properties 
of the reduced characteristic equation are investigated. For the reason 
of simplicity and convenience most of the examples in this chapter take the 
coefficients themselves as variables. For the case of the adjustable para- 
meters scattered in all the coefficients, the analyses are the same. In 
fact, where the analysis of the reduced characteristic equation alone is 
concerned, it does not make any difference about what type of compensator 
is used. The analyses in this chapter also assumed that no restrictions 
are imposed in the variables, consequently no restrictions are imposed on 
the arbitrary roots. 

2-2. Coefficients as variables. 

Treating some of the coefficients of the controlled characteristic 
equation as variables provides a convenient and relatively systematic 
means for the purposes of analyses. However, the actual applications 
have some difficulties. Theoretically, the variables can be transformed 
analytically or graphically from one set of variables to another by the 
root-coefficient relation and the definition of variables. Thus the 
variables which originally are coefficients of a characteristic equation 
can be transformed to another set of variables such as the actual adjust- 
able parameters. However, the transformation may be very complicated 
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because of the inherent nature of the algebraic equations. In actual 



design, it is unnecessary, because the actually adjustable parameters 
can be treated as variables directly by the partitioning processes 
described in the last chapter. 

Mitrovic's method treats two coefficients as variables. The applica- 
tion of Mitrovic's method to the case in which the actually adjustable 
parameters enter some of the presumed fixed coefficients requires that 

transformations or approximations must be made either graphically or 

2 4 5. 

analytically. Some of the work to this end has been done. 9 9 
2-3. Effectiveness of the adjustable coefficients. 

The effectiveness of stabilizing ability of adjustable coefficients 
are discussed in this section. Since only a qualitative analyses, i.e,, 
the tendency of stabilization is required, the root-locus method is 
applied to the reduced characteristic equation. 

Consider a fourth order characteristic equation (2-1' the effectiveness 



of adjustable coefficients are analyzed in several cases t 

Case (1): adjustable. Since only I $2 is adjustable, then the 



5*4 6j5"4- 4 B,S 4- 8 0 ~0 



(2-1) 



roots are constrained by B, and B • Assume the reduced characteris- 

J 3 1 o 

tic equation is 



S 5 4 C 2 5 z 4 C,5 + C 0 — O 



( 2 - 2 ) 



and the arbitrarily chosen root is -S^ 9 then 

®3 = C x + Si 



B 



1 



C 0 4 C.S 
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From the above equations, solve for C^, and C qS> obtain 

Q = B, - 5, 



(2-3-1) 



c, = y, ( b,-c„) = 




( 2 - 3 - 2 ) 

(2-3-3) 



Substitute (2-3) into equation (2-2), one obtains 

S 5 + ( B,- S.ji" + (-37) 5 - -§= C S - S,) = O (2-4) 

Equation (2-4) is the reduced characteristic equation with the coeffici- 
ents expressed in terms of the arbitrary root -S. , Assume S. has been 

1 1 

chosen, the tendency of the root locus of equation (2-4) by varying B q 
is to be investigated. Equation (2-4) in root- locus form is 




(S-S.) 



s[sN ( b 3 





(2-4) 



Consider equation (2-4) and assume then there is a zero in the 

right-half plane and the angle-criteria of the root-locus is 2.7LH (n = 0, 
1, 2,) instead of m TL (m = 1, 3, 5...). The general root-locus pattern 
is shown in Fig. 2-1, in which and R£ are the roots of the equation 




Fig. 2-1 Root-locus of equation (2-4) Case of 
1*2 adjustable. 
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Case (2): adjustable* Now consider the same fourth order 

characteristic equation (2-1) in which only the coefficient of the 3rd de 
rivative term is adjustable* By the same manipulation as Case (1), the 
coefficients of the reduced characteristic equations are obtained as 
follows : 



c, =. 6 ,/s , - B c /sf 



(2=5-1) 



— Bj,/ S, - B,/5, 2 -t B 0 /s * 



( 2 - 5 - 2 ) 



C 0 = Bo/ S, (2-5-3) 

Substitute (2-5) into (2-2), the reduced characteristic equation be- 



comes 



s > +(-lr;-ji)s 1 +( J si)s +-js(s‘-s,s +s>J=o 

In root-locus form for B q as variable s (2-6) becomes 



( 2 - 6 ) 



3; 



^ ( 5 ~ 5,5 + 5, J 



s[s‘ + (f -fns + f] 






(2-7) 



Assume — \ - ■ 0, equation (2-7), has two complex zeros in the 

right half plane and the angle-criteria of the root-locus is + 180°. 
The root locus pattern is shown in Fig. 2-2. 




o 






CV 




Fig. 2-2. Root-locus of equation (2-7), adjustable 
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Case (3): and adjustable. 

Consider the same fourth order equation but with and adjustable. 
Since there are two coefficients adjustable, the reduced characteristic 
equation is of order 2 # Follow the same derivation, the coefficients of 
the reduced characteristic equation are: 



C,= -sTTsJ B,-5A-c 0 ) 

Bo 

5, 

where -S^ and -S 2 are arbitrarily chosen roots and 
are assumed. The reduced characteristic equation is 



( 2 - 8 - 1 ) 



( 2 - 8 - 2 ) 



5 + ^( 6 ,-^ 




(2-9) 



In root-locus form for B as variable is 

o 



V 



(5,+ 5i) J _ 






(2-10) 



This equation again has one zero in the right half plane and the angle- 
criteria of the root-locus is lILf 1 (n = 0 8 1, 2„...). The general 
root-locus pattern is shown in Fig. 2-3. 




Fig. 2-3. Root-locus of equation (2-10). and 
adjustable. 
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Case (4): and adjustable. 

In this case, the reduced characteristic equation is of order 2 
and the coefficients from Table 1-2 are as follows: 



C, = T 



B. 



3 , 5 , 



(s,5 x y 



a 0 



( 2 - 



c< 



Bo 

3,5* 



( 2 - 



In root- locus form* the reduced characteristic equation becomes 



( 5 | ) o t s — \ 

D oW 5 ,^) 



3,5, 



~ / 



( 2 ° 



Equation (2-12) has one zero in the right half plane and the angle- 
criteria of the root- locus is ZKf\ (n = 0 S 1 9 2 B ...). The general 

root-locus pattern is shown in Fig. 2-4 when 0. 




Fig. 2-4. Root-locus of equation (2-12). 1*2 and 

adjustable . 



11 - 1 ) 

11-2) 

12 ) 
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Case (5): and adjustable * From Table 1-1. The coefficient 

of the reduced characteristic equations are as follows: 

C, = B 3 -(5 J + 5 i ) (2-13-1) 




(2-13-2) 



The reduced characteristic equation in root-locus form is 






Bo 



3 [ 5 + ( dj-5, -4)J 



= -I 



(2-14) 



The general root-locus is shown in Fig. 2-5. 




Fig. 2-5 Root-locus of equation (2-14). 
and B 2 adjustable. 

Case (6): adjustable. In the same way 9 the reduced characteristic 

equation in root- locus form is obtained as follows: 



Bo 

5/ 

s[s x +(Bj-S,)S + 5jS,) 



(2-15) 



The general root-locus pattern is shown in Fig. 2-6. 
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Fig. 2-6. Root-locus of equation (2-15) 
adjustable. 

Case (7): Bp and adjustable- 

In this case three roots are arbitrary,, only one root is constrained. 

Assume S, is the constrained root, and define C = S, , then 
4 o 4* 



3 0 — S.S.SjCc 

C„= Se/s.iiil 



The reduced characteristic equation is 

Bo 

5 + 



o 



If no restrictions are imposed on the adjustable coefficients, then for 

any choice of Sp and 3^ 0, ^ 0, 0), the system can be 

stabilized for any value of B * 
y o 

From the above analyses and the general pattern of the root- locus 
three conclusions can be drawn- 

(1) If the adjustable coefficients are not in the lowest order 
or not in sequence, as the cases (1), (2), (3) and (4), right half plane 
zeros are introduced in the root-locus plot. Consequently, the stabiliz- 
ing ability is less than that if no right half-plane zeros would be 
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introduced. 



(2) If the adjustable coefficients are in the lowest order 
and in sequence,, as the cases (5), (6) and (7), no right half-plane 
zeros are introduced. Consequently the stabilization is more effective, 

(3) The most effective scheme is the case where all of the 
coefficients are adjustable as case (7)* The most Ineffective scheme 

is the case where only the highest coefficient is adjustable' as case (1), 

In general s if there are 9C r sv unadjustable coefficients which have an 
order lower than the highest order adjustable coefficients, then the 
reduced characteristic equation in root-locus plot for as variable 
has * # r n zeros in the right half plane. Because of this property, it seems 
that if the adjustable coefficients of a characteristic equation are not 
the lower order and not in sequence, the stabilizing tendency of these sets 
is less than for the case of adjustable coefficients which are of lower 
order and in sequence. Consider case (5), in which and B, t are adjust- 
able, the reduced characteristic equation is of the second order and there 
are no zeros in the right half plane, consequently the system can always 
be stabilized for any value of K. However, for cases (1) to (4), it 
cannot always stabilize the system. In general, for any order system 
if the reduced characteristic equation is of second order and all adjust- 
able coefficients are of lower order and in sequences, then .the system can 
be stabilized for any forward gain K. However, this is not the case if 
the adjustable coefficients are not of lower order and in sequence. More- 
over, the more zeros introduced in the right half plane in the process of 
derivation of the reduced characteristic equation, the less the tendency 
to stabilize the system. This implies that the more the unadjustable co- 
efficients which are of lower order than the highest order adjustable 
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coefficient, the less the effectiveness of the control parameters. This 
is the inherent property of algebraic equations. It is because of this 
property that many approximations can be made by ignoring those control 
parameters which affect only the coefficients of higher derivatives. If 
the design of a control system is based on the characterist ic equation, 
those coefficients of lower order derivatives should be considered first 
if it is possible. 

This analyses also explains the well known fact that a cascade filter 
is more effective than lower order derivative feedback, especially for 
high order systems. For a cascade filter, the filter parameters are 
scattered in all coefficients of the characteristic equation in the nature 
of algebraic equations, while derivative feedback enters some of the co- 
efficients only. This analyses also explains the fact that tachometer 
feedback is good for low order systems, but may even hurt the performance 
of high order systems. 

2-4. Root-Region of Stability. 

The reduced characteristic equation is ,a function of the arbitrary 
roots. Mathematically for any choice of the arbitrary roots, there ex- 
ists solution of the equation. However, in control systems, the roots 
of the characteristic equations are required to be in the left half plane, 
this requirement therefore confines the arbitrary roots in a certain region 
in the left-half plane. This region here is called "region of stability' 1 
and may serve as a guide for the choice of the arbitrary roots. 

If the number of variables is one, this region is a line or lines 
in the S-plane (root-locus). When the number of variables is more than one, 
this region is an area or several areas on S-plane. The boundary of this 
region can be calculated by Routh*s criterion subject to the reduced 
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che racteristlc equation. The pattern of 4th and 5th order equations for 



h - fi j h ~ f j> 

and Fig. 2-8. 



and -f- f - as variables are shown in Fig. 2-7 
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(b) r,,.f 3 
variable 




Fig. 2-8 Root region of a 5th order equation 



(c) f,,f 3 
variable 
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In Fig- 2-7 and Fig- 2-8* if a pair of complex roots are in the 
stable region* then all the roots are in that region- If a pair of 
complex roots are in the unstable region* then there is at least one root 
in the right half plane- If a pair of complex roots are on the stability 
boundary line* then there are at least one pair of complex roots on the 
j 6 0 axis or one root at the origin. The arbitrarily chosen roots* 
therefore* must be confined in the stable region- The calculation of 
the boundary is illustrated in the following example. 

Example 2-1; Given a fifth order equation 



Case 1; Assume and are variables* then the controlled 

characteristic equation becomes 



5 r 4 /(cS^-h US S* + £2.os Z -h 1300 5 4 2006 =0 ( 2 - 16 ) 



5 * + I kS 4 4 f 3 S* 4 fiS 2 4 ISoos + Zooo — O 



(2-17) 



From Table 1-2 ; 



c, = 8+- (S.+5J 



(2-18-1 ) 




(2-18-2) 



B 0 



C o S,5 Z 



(2-18-3) 



where -Sj and -S^ are the arbitrarily chosen roots. Assume these roots 

are confined on a constant ^ line, then 



S 2 — j — CO n 
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Substitute into (2-18) 



c 2 = -2fUJ n 



(2-19-1) 



1 r . if ] 

C. = <** l 8,- 



(2-19-2) 



Co — 



J3° 

<*>n 



(2-19-3) 



where = 16, = 2000, B^ = 1300 from the given equation (2-17) . 

The Routh's criterion for the 3rd order reduced characteristic 



equation is 



C,Q > Co 



(2-20) 



Substitute (2-19) and the values of and B^ into (2-20), 

and manipulate, the Routh's criterion in equality becomes 



Z y (o f^n 



iOn 



< i d.8 + af 



(2-21) 



Let the left side of the (2-12) be -j- , 

f = + ^ 

The -j- versus Uj n plot for given value of 



then 



f is shown in Fig. 2-9. 
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From the inequality (2-21), if oJ, £: 0U n — , then all the roots 

of the reduced characteristic equation are on the left -half plane. CO / 
and define the stability boundary of the arbitrarily chosen roots and 
are the solution of the following equation 



Z.(o 4 = 13.8 + 8 ^ 



( 2 - 22 ) 



Solutions of (2-22) for O J ^ I are tabulated in Table 2-1, and 

the boundary for root region are plotted in Fig. 2-10. 



Table 2-1. Solutions of equation (2-22) 



f 


CO, 


CO± 


1 


3.7 


6.62 


0.9 


3.28 


7.53 


0.7 


2.4 


10.5 


0.6 


2.05 


12.9 


0.5 


1.7 


14.3 


0.4 


1.35 


17.9 


0.3 


1 ,05 


23.95 


0.2 


0.65 


36.1 


0.1 


0.35 


72 



Case 2: f, f variable 

' i y i 

Since the coefficients of the reduced characteristic equation are 
expressed by the arbitrarily chosen roots, the mapping contour can be 
chosen arbitrarily. In order to demonstrate this, let the two arbi- 
trarily chosen roots be on a constant Online, then 



5 , + 5 * = (2-23-1) 

= <9^ -t LO^ (2-23-2) 
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joJ 




Fig. 2-10 Stable root region of a 5th order equation with 
variable of example 2-1. 
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Substitute (2-23) into the coefficients of the reduced characteristic 
equation in Table 1-2, one obtains 

C z = It, - 2.0^ (2-24-1) 



c, == A-2 0 -( - (of + o/j (2-24-2) 

C 0 = ZOOCy^C + tO*J (2-24-3) 

The Routh's criterion is C Z C, 3 Co Substitute (2=24) into the 
above inequality and manipulate, obtain 

(|k- 2 <v) |28 -(of iW jJ (<?c f OjJ Z.OOO (2-25) 



The boundary value of CO for the inequality (2-25) are solved by the 



same procedure as Case 1 and tabulated in Table 2-2. The stable root region 



is plotted in Fig. 2-11. 



Table 2-2 Stable root region of a 5th order equation 
with f, and variable. 



(K 


Stable 


region of 


lO 


0 


0.975 


< 


cu 


< 


11 


.29 


2 


-8.57 




uJ 


< 


+8 


.57 


4 


-6.62 


c 


uJ 


< 


+6 


.62 


6 


-5.9 


< 


uJ 


< 


+5 


.9 


7 


-6.3 


< 


UJ 




+6 


.3 


7.5 


-5.9 




uJ 


< 


+6 


.9 



-0.975 > o 0 > - H.29 



Case 3: -f, f, variable. 

1 I J 1 J 

In this case, from Table 1-3, the coefficients of the reduced character- 



istic equation are: 

Q — 3 4 -(3, + 5 z ) 



(2-26-1) 
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Fig. 



unstable 

region 




Stable root-region of a 5th order 
equation with -f, and f * variable 
in example 2-1. 



-f- (Tv- 
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5 , 4 =*p O 



( 2 - 26- 2 i 



C,= ijs, [B i -C.-C,SA] 

Co = Bo/s.Si 
Assume Sj and S 2 are on a constant j* line. Then 

s , + 5* = 

S,5 Z = U)n 



(2-26-3) 



(2-27-1) 

(2-27-2) 



Substitute (2-26) and (2-27) into the Routh criterion inequality and 
manipulate, obtain 



-/6cjJ +ij , cq i > )uJ n < £000 < 2 - 28 > 

This inequality is a function of Uj n U P to the 5th order, the solution 
is more easily found by a graphic method. The function of the left side 
of (2-28) for J as a parameter are plotted as a function of C <J n > 

boundary values of C<J n are evaluated as shown in Table 2-3. The root- 
region of stability is plotted in Fig. 2-12. 



s 


S^ab 1 


e region 


0.1 


2.2 


< 


U) n <5.45 


0.3 


2.3 


< 


U) n < 6 


0.5 


2.4 


< 


tU n < 7.1 


0.7 


2.5 


< 


LO n <10. 7 


0.9 


2.6 


< 


UJ n < 8.6 


1.0 


2.68 < 


OJ n <7.3 



Table 2-3. Stable region of 5th equation with 
, .j: variable. 
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Figure 2-12. Stable root-region of the 5th order equation with 
j. , -j variable in Example 2-1. 
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In equation (2-26-2), 3 was assumed. This implies 

j 0. For ^ as 0, the coefficient of the reduced characteristic 
equation is undefined. Rearrange equation (2-26-2) 



C,( 5, + 3 2 J — Bj~ 
For ^ = 0, 5, + 5^ = 0 then 

B z - C 0 — uj l — O 
B 4 ^- B.ujS 0 o =0 



O ~ C, 5,3^ 



( 2 - 2 S'J 



consider equation (2-285, the coefficients and are unadjustable 

coefficients of the given characteristic equation, therefore, the solu- 
tion of (2-28) defines four points ijC 0, » and 4 on t * ie 

axis instead of a region. Those four points are convergent points for the 
stability boundary, i.e., if there are any roots on the j 00 axis for this 
type of control, the roots must be either 1 , or + j uJ^ regardless 

of where the other roots are. Consider Fig. 2-13 which is the same as 
Fig. 2-12, but repeated for convenience. By the definition of stability 
boundary, if there are any roots on the boundary, there are at least one 
pair of roots on the Jo 0 axis or a root at the origin. Assume a point 
P (inside the stable region or the unstable region) is a root of the equa- 
tion, as P approaches the boundary in any manner whatever, equation (2~28) r 
implies there are either a pair of complex roots which approach the points 
1 j 6J, or » or a r °ot which approaches the origin. In case one 

and case two of this example, the roots on the jOJ axis have an interval, 
and the location of the roots in this interval depends upon where the other 
roots are. In this case, the roots on the j OJ axis are independent of the 
other roots. In general, if the variable coefficients are not in sequence, 
then the roots on the j (jJ axis have this property. In short, for this 
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Fig. 2-13 Convergence on the j Cd axis 
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type of control if the variable coefficients are chosen such that 

there are roots on the >6 O axis, then those roots on the J CO axis are 

completely determined by the fixed coefficients of the given equation., 

/ 

In this example, substitute B^, and into equation (2-28), 
one obtains 



l O, = 2.11 j S.Z 

The discussions above are for complex roots. For real roots the 
analysis can be carried out just by the same procedure. 

2-5. Dominant Root Region: No Zeros Close to the Origin, 

Since most of design is based upon a pair of dominant complex roots, 
a procedure of evaluating the "dominant root region" is to be discussed. 
The definition of dominant roots depends upon the degree of dominance and 
other requirements. It can be defined as one wants. In case of no zeros 
close to the origin and if the pair of dominant complex roots are so de- 
fined such that the negative real part of all other roots (complex and 
real) of the original characteristic equation must be greater than that 
of the presumably dominant roots, then the Routh’s criterion can be 
applied to the reduced characteristic equation by a conformal mapping 




i<A 



(a) S-plane 

Fig. 2-14 Conformal mapping of S = p-a 
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Consider Fig. 2- 14(a), assume the pair of arbitrary complex roots 
(- S ^ ) and ( - S 2 ) are the dominant roots as defined in the last paragraph 
and also assume all the other roots of the original characteristic equa- 
tion must have the real parts which are greater than ,J a 98 (a is positive 
real number and greater than or equal to )» then all the roots of 

the reduced characteristic equation must be to the left of the line S - -a 
in the S-plane* By the linear transformation S - p-a, the region to the 
left of the line S = -a in the S-plane is transformed to the left -half plane 
in the P-plane. Therefore the Routh e s criterion applied to the p-plane, 
implies all the roots of the reduced characteristic equation are to the 
left of the line S = -a in the S-plane. In addition, to put all other 
arbitrarily chosen roots (there may be more than two arbitrarily chosen 
roots) also to the left of the line S = -a, then all the roots of the 
original characteristic equation except a pair of complex dominant roots 
are to the left of the line S * -a. Here the line S ~ -a was assumed to 
be fixed if the dominant roots are fixed. In order to find the dominant 
root-region which satisfy the definition, °a° must be a variable as the 
dominant roots are varying. Therefore the application of Routh e s criter- 
ion to the p-plane defines a region in the S-plane if such a region exists. 

To obtain the transformation from S-plane to p-plane, the process is 
just a substitution of S = p-a into the reduced characteristic equation. 
Consider a third order reduced characteristic equation 

S 3 + C z 5 2 + C,S 4 C c = O (2-29) 

substitute S - p-a into the above equation, one obtains 

C 2 (P-a/-f- C.(p-Q) +C Q = O 
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manipulating, it becomes 



U z -3a.) p z + ( C, - zac z +3 clJ p +(C 0 -C,aiQa*-a!J (2 ' 30) 

=rO 

Equation (2-30) is the transformed reduced characteristic equation in the 
p-plane. Now define: 

= C Z ~3C L (2-31-1) 

-P, =■ C,- 1 <\C L 4 3 a . 2 (2-31-2) 

D a = C„ - C,a. + Q of - a. 3 (2-31-3) 

Where D’s are the coefficients of the reduced characteristic equation in 
p-plane, then equation (2-30) becomes 

5* + P i S i + X>,s + J> 0 = O (2-32) 

The Routh criterion of equation (2-32) is 

D, D 0 (2-33) 

As the coefficients D’s are functions of the arbitrarily chosen roots and 
the definition of ’a 1 , the inequality from Routh 1 s criterion defines a 
region on the S-plane directly. 

The coefficients D's for a nth order reduced characteristic equation 
and its Routh' s criterion inequality have been formulated and tabulated 
in Table 2-4. Since the reduced characteristic equation of any order 
system and any number of variables can be derived,, the process of evalua- 
ting the dominant root-region for the definition described above can be 
applied to any order system and any number of variables. 
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Table 2-4. Transformed coefficients D's of the reduced characteristic 
equation "a" is the mapping contour, C # s are the coefficient in the S-plane 



order of the 
reduced 
equation 


transf orraedcoef f . 's 


Rough's 

criterion 




Px = 1 
^1 Si Cj 

P* = C;-CA 4 <X 


J>,>0 

J>o >° 


3 rd 


p, = 1 

Pi. — C z - 3 ^ 

Pi = c, ~-zoCl- 15 cZ 
p 0 = c 0 -c,a 4 ^-0.3 


J>,>0 

P^>0 

M >A 




P4 ~ 1 

P3 = C } - 40 - 

Di = C 2 - 3 dCj iba? 

P, — c, - 20 Cj fSoCj-40! 

P 0 = Co - ac, i <x l C 1 - clCj t of 






\ 

n 


J>„ = Cn 

p„., = c„-, - n«.c„ 

Pm= C^-Cn-OaC^+^V^ 

P„_ 3 = -Cn-2;aC 0 . 2 . + L; ^r' 

J _ n(n-i)in-z) _} ^ 

h 

1 - 2-3 

-Ri-t - C n . + -c n-3) a C 0 - 3 1 C 

- l Cn t 

\ • 1 •>* + 

# < 

• 

4 


:n 

Crv 

D-2 


t 



I 
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Fig. 2-15 Plot of inequality 2-3/. 
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Example 2-2. Consider the same characteristic equation given in 

Example 2-1 9 namely - 



5 > + -f 6jS* + t,S 4 6 C — 



(2-34) 



where B, = 16, B_ = 128 s and B » 2000* The coefficients of the reduced 
4 3 o 

characteristic equations are: 



2. 

'n 



(2-35-1) 

(2-35-2) 

(2-35-3) 



C ^ — / (o — 2 j ^ n 
C | — U8 - 3ZU)n + (4f- 'J 
C 0 — ZOQo/ UJ* 

where £ and uU r define the dominant roots^ since there are only two 
arbitrarily chosen roots. Assume a » 9 i.e.j all the other roots must 

have real parts greater than or equal to ^CaJ n 9 Substitute equation (2-35) 
into (2-31) the transformed coefficients of the reduced characteristic 

equation become: 

p z =. Ib-SrflOn (2-36-1) 

P, =■ I 2 8 - fctf UJ„ + ( >>f- 0 idn (2-36-2) 

J) o — 2000/f^ - I2& fU) n +48f Z C^ 1 i - j 5 ( Tf - 0 uJ,f (2-36-3) 

By substitution of (2-36) into the Routh’s criterion DjDj > D q and manipul- 
ating, one obtains: 

l&LfC0 n -U> (23f-i>J Z - *Y$ ( I- 12 f) ^ 2o<^6 ^ 2 “ 37 > 

Let the left side be -j- ( f col , choose ^ and plot -f- versus CU n as shown 

in Fig. 2-15, the dominant roots (defined by a = f{J t ) region which satisfy 

the definition for £ = 0.3, 0.5, 0.7 are evaluated from Fig. 2-15 as 

follows : 
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6.3 



f - 0.3 1.4 ± <\) n e 

f =0.5 1.65 UJ n £: 5 

j 5 = 0.7 1.9 £ u) n ^ 4.8 

This region is plotted in the S-plane as shown in Fig, 2~16 s in which 
the stability boundary is also plotted. 

In the dominant root- region, assume the dominant roots are chosen to 
be -S l8 -S 2 =1.89 + j 6; ( f = 0.3, t^,«6.3) 

the secondary roots evaluated from the reduced characteristic equation 
are : 

-S 3 , -S 4 = -2.8 + j 1.0, ( £= 0.94 9 ^ n ~ 2.82) 

-S_ e -6.6 
5 

It shows that the secondary roots have real parts which are greater 

than that of the dominant roots. 

If the dominant roots are chosen as 

-S v -S 2 = 2.5 + j 4.32 s ( f » 0.5, W, = 5) 

the secondary roots are- 

-S 3 * -S 4 = 3.56 + j 2.76, ( f « 0.79, t J n m 4.5) 

-S e = -3.88 

It shows again the real parts of the secondary roots are greater than 

that of the dominant roots. 

In the above two choices of the dominant roots , all the secondary 
roots also have greater damping ratio ^ . However the dominant root region 
evaluated by this procedure doesn°t imply this restriction to the secondary 
roots* but this restriction can be obtained by choosing large value of "a" 
relative to by inspection of the shape of the region. 

In Fig. 2-16, it can be seen that the dominant root-region is inside 
the stability root-region. The secondary roots are not necessarily inside 
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the stability boundary. The max. damping of the dominant roots can 
also be evaluated from the dominant root -region,, 



Example 2-3. Given a 6th order characteristic equation 

s‘°+s ? s i '+ 8^“ + f s s J + ts'-tf.s ■+ 6„=o 

Q 

where f { , f „ are adjustable, - 184 y - 10840 and = 3x10 . 

Here three coefficients are adjustable, three roots can be chosen 
arbitrarily. Assume the requirements of the dominant roots is such that 
all the secondary roots must have the real negative parts at least twice 
that of the dominant roots, then the mapping contour is a = - In 

order to meet the requirements of the secondary root, the third arbitrarily 
chosen root must be set at least to the min. value, namely ^ 2 f # 
Assume - 2^t0^ then the three arbitrarily chosen roots are: 



"3, j ~ J f ~ " dominant 

- s 3 — 

From Table 1-5, the coefficients of the reduced characteristic equation 
are : 



C z = B s -( + 

C, = %- + +S.S.J 

C„ = 6 ,/ASj 

Substitute S^, S 2 and into (2-37) and manipulate, obtain 



(2-37-1) 

(2-37-2) 

(2-37-3) 



c a -- 3^-4-f^n 

C , = - 4? Bj.U u + ( uf- 0 UJn 

C 0 = B„/2?U> 0 5 



(2-38-1) 

(2-38-2) 

(2-38-3) 
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Fig. 2-16. Dominant root region of a 5th order equation with t,, fi variable. 
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fig. 1/ 



Dominant kooc Regions ot Example z-j 
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Substitute (2-38) into the coefficients of the transformed reduced character- 
istic equation in Table 2-4, obtain 



£> 2 “ d s ~ 10 f -A, 

P, : - 8f i + (lof- I ) l0,~ 

V — Bo/zfuJn - + > 2 f Z B s UJ* - 2J (z4J - >J u)n 

The Routh*s criterion D 9 D > after manipulation becomes 

z l o 

Bj + 5JUJ n + B s ( 103 f- ') UJ* - Sf ( 4+f - l) j > 3/1 f 

Substitute the numerical values of B es B, and B , it becomes 

5*4 Q» 



> ! (r . , Z \ r 2. \ 3 V hr* to* 

: ISLiiO ^ 4 -Oil),) | (2*39) 



For J =0.3, 0.5, 0.7 and 0.9, the region of which satisfy the in 

equality (2-39) are evaluated and plotted as shown in Fig. 2-17. 



f 


= 0.3 




7.5 - 30 


f 


= 0.5 




7 ‘= <^-19.5 


f 


- 0.7 




6.8 ^ tj„f=13 


f 


= 0.9 




None 


Within the dominant 


root region, 


if the 


dominant roots are chosen as 


_S 1> -S 2 = -9.5 + 


J 16.4° 


< t 


* 0.5, U) n = 19) 


then: -S 3 = -2 ( 


» -19 







The other three roots evaluated from the reduced characteristic 
equation are: 



-S 4 , -S 5 = -20 + j 25 ; ( = 0.97, = 20.6) 

-S 6 « “106 

If the dominant roots are chosen as 

-S 1 , -S 2 = -8 + jl3.8, ( f = 0.5, cJ n = 16) 

The secondary roots are 
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The evaluation of the dominant root region is not necessarily carried 
out point by point, only those points of interest are evaluated. And 
furthermore, only the boundary points which satisfy the Routh’s criterion 
are of significance. The Routh's criterion formed in this process is a 
polynomial of LO n in any case, therefore for higher order reduced 
characteristic equations there may be more than one region which satisfies 
the criteria. But only the region which is closest to the joJ axis is of 
interest, because the roots are assumed to be dominant. 

2-6. Modified Dominant Root Region - One Zero Close to the Origin. 

The assumption that the response of a system can be characterized by 
a pair of complex dominant roots is based on the approximations that all 
other roots are either too far to the presumable dominant roots or close 
to a zero. Consider Fig. 2-18, the roots (»S^) and are said to be 

dominant 



jcJ 




'ii 



X 




- 5 , 



+<K 



- S, 



X 



Fig. 2-18 General Closed Loop Pole-zero Configuration. 
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because the other roots (-S,-)* (-S^) and (~S^) are far to the left 
and (-S^) Is close to a zero (-Z^) so that they contribute negligibly to 
the transient response* If no zeros lie to the right of the dominant 
roots (close to the origin) * then all the secondary roots are required to 
be far to the left if the presumed dominant roots are dominant * If there 
is a zero close to the origin* then one root to the right of the domin- 
ant roots has negligible effect on the dominance if that root is close enough 
to the zero. Therefore* the dominant-root region defined in the last sec- 
tion must be modified when there are zeros close to the origin* Usually* 
the zeros of a system are known. If the known zero is far to the left of 
the dominant roots* then the dominant root region defined in the last sec- 
tion is still valid* When the zero is known to be close to the origin* 
the modification of the dominant root region is required* 

Consider Fig* 2-19(a) which is a third order system with a cascade 
compensator* Assume the compensator is a lag network* Let (~sp and 
(- S 2 ) be the arbitrary roots (also the dominant roots in this case)* then 
(Z) and (p) are functions of and (Equation 1-22)* The general pattern 
of the root locus is shown in Fig* 2-19(b)„ Since and $2 are arbitrary 
roots* then from the tendency of the root locus* the one only root which is 
possible close to the origin is (-S^). If the arbitrary roots are chosen 
within the dominant root region (point B in Fig* 2-19E) defined in the last 
section and the pole (p) and zero (Z) of the compensator are adjusted ac- 
cordingly* then (- S ^ ) and (-S^) are to the left of the dominant roots (Fig* 
2-lSb). If and $2 are chosen on the boundary of the dominant root re- 
gion (Point C in Fig* 2-19E), then the root (-S^) is forced to be in a 
position as shown in Fig* 2-19C. As the arbitrary roots move toward the 
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Fig. 2-19a Cascade filter. 
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Fig. 2-20 Dominant root region (a = and Modified Dominant 

Root Region 
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origin and so does the zero. When the arbitrary roots reach the 
boundary of stability (Point D in Fig. 2-19E), (-S^) reaches the origin 
and the zero and pole of the compensator also reach the origin as shown 
in Fig. 2-19D. From this analysis, it can be visualized that if the 
dominant roots are chosen close to the stability boundary and within 
the stable region, this implies a dipole close to the origin. Therefore, 
a region which is close to the stability boundary and between the boundaries 
of stability and dominance (which are defined in the last section for a = 
^U)n ) is also a dominant root region. This region is given the termin- 
ology n modified dominant root region** in order to distinguish it from the 
definition of the last section. These two dominant regions are shown in 
Fig. 2-20. 

The pole-zero configuration discussed above is often the case for 
systems of high velocity error constant (K^) . For very high velocity 
error constant, a dipole near the origin is necessary. 

The analysis above is based on the tendency of the root-locus and 
the known zero. The same analysis can be applied to other cases. 

In this case the dominant roots should be chosen either close to the 
stability boundary or within the dominant region (a = ). Choosing the 

dominant roots in the vicinity of the dominant boundary (a = j Oj ^ ) 

should be avoided when a zero is close to the origin because the root is 
not close enough to the zero and it may cause excessive overshoot. If a 
zero is known to the right of the dominant roots, a root close to it is 
necessary. 

Now consider the pole-zero configuration of Fig. 2-19D which is the 
situation of stable limit. In this condition, the pole and zero of the 
compensator are at the origin, therefore the presumed dominant roots -S^ 
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and -S^ are on the original uncompensated root-locus (K as variable). 

This provides another means to evaluate the stable boundary, namely the 
stable boundary is the uncompensated root-locus itself in this case (lag 
compensator). This boundary line is one of the dividing lines of the 
root relocation zone discussed in reference 7. For high order systems 
there are several branches of the stable boundary line, but for the pur- 
pose of dominant region only the one which is closest to the jCO axis is 
of interest. From the root-locus it can be visualized easily that the 
dominant roots are within the region which is the closest to the jUJaxis. 
From this analysis, it concludes that for lag compensator the region 
closest to the nearest branch of the uncompensated root-locus is the 
dominant root region. And the dominant boundary of the definition a = fCO n 
provides a guide of how close they should be. 

Example 2-4. In Fig. 2-21, is not to be reduced. Find the dom- 
inant and modified dominant root region: 



o 



*(3+£) 




(O 


(s + pj 




Sif+ZS-H.Z) 



Fig. 2-21 Example 2-4. 

The expressions for the dependent variables Z, p, C q and have 
been derived in (1-22) and (1-23). Substitute » ^ 1^2 = 

(jJ^ , B q = 6, = 1.2, B 2 = 2 into (1-22) and (1-23), one obtains 
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C 0 = [(1-4. fj<-u + +jH] 6 / ( ^UJ*-^50J S n -(o) 



(2-40-1) 



P = ^‘c»/fc 



(2-40-2) 



C, — Z + p ~2fOO n 



(2-40-3) 



CX = [ td n C 0 ~ l-Zp] / (o (2-40-4) 

From Table 2-4: 

P, = c,- aa 
J>o= Q-C.a+a 1 

For a = 

13 1 (z-4-l-l) 

Po = c„- H,c, + *'< < 2 - 41 - 2 > 

The Routh's criterion for D- and D is D, .> 0 and D ^ 0. 

1 o 1 o 

The stability boundary line (root-locus of the uncompensated system) 
is shown in Fig. 2-22a. In order to evaluate the a = ^ UJ n dominant 
boundary line, C q and must be calculated. For ^ = 0.3 and ^ = 0.5, the 

results of computation from equations (2-40) are tabulated in Table 2-5a 



and 


Table 2-5b respectively. 


The dominant 


boundary 


line for a 




is 






Table 2- 


5a. Example 2- 


-4 for f = 


0.3 








Co 


c, 


P 


©< 


i 




A 


0.4 


0.65 


1.777 


0.01735 


0.0698 


0.0249 


+ 


+ 


0.5 


0.473 


1.72 


0.0197 


0.0914 


0.216 


+ 


+ 


0.6 


0.278 


1.657 


0.01662 


0.118 


0.141 


+ 


+ 


0.7 


0.0539 


1.584 


0.0044 


0.1322 


0.0335 


+ 


- 


0.8 


-0.085 
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Table 2-5b Example 2-4 for ^ = 0.5 



U) n 


Co 


c, 


P 


04 


Z 




K 


0.2 


0.81 


1.8 


0.0054 


0.0543 


0.0995 


+ 


+ 


0.3 


0.615 


1.7 


0.00923 


0.0545 


0 1695 


+ 


+ 


0.4 


0.418 


1.611 


0.01115 


0.0685 


0.1625 


+ 


+ 


0.5 

0.6 


0.213 

0 


1.509 


0.00886 


0.079 


0.112 


+ 


' 


Shown in 


Fig. 2-22a. If the 


dominant roots 


f = 0. 5 S 


U)= 0.55 


are 


chosen 


which is 


close to 


the stability boundary , 


the closed 


loop pole 


-configure- 


tion is 


shown in 


Fig. 2-22b. 


For this choice of the 


dominant 


roots 


, the 



compensator is 




o.o Sf Q^oU) 
S + cooS^-3 
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CHAPTER III 



FUNCTION OF VARIABLE PARAMETERS 



3-1 Realization Problem: 

The result of the partition process employed in this method is to 
separate the dependent variables into two parts: the coefficient of the 

reduced characteristic equation and the variable parameters (parameters 
of the compensator). The design is carried out essentially by two in- 
dependent procedures. One is to choose appropriate closed loop pole- 
zero configuration from the reduced characteristic equation and the other 
is to realize the compensator according the chosen pole-zero configuration. 
Consider a cascade compensator of the following transfer function. 



If there are no restrictions to the type of compensator, the parameters 



Their solutions from the partition functions are always real. For the 
function of (3-1) with all real coefficients always can be realized either 
by a passive network or an active network. However, if the poles of the 
compensator are required to be on the negative real axis, then (3-1) be- 
comes 




K c ( 5*4 Q-,5 f (Ko) 
C sN b,5 4 b 0 ) 



(3-1) 



a^> a^, b^ and b^ are linearly dependent in the root-coefficient relations. 




(S-i-p.XS + p*) 



(3-2) 



From (3-2) and (3-1) 



Pi + P 2 . = b, , 



R R = bo 



(3-3) 



In equation (3-3), and P ^ are quadratic functions* For certain choices 
of the roots the solution of P^ and P^ may be complex numbers* although 
and jo o are always real* Therefore the choice of the arbitrary roots 
are not only necessary in the dominant root region but also necessary in 
the region which gives a physically realizable compensator. For certain 
kinds of specifications and restrictions* the compensator may be forced 

g 

to be an active network * but for most of the systems* the designer should 
try a passive network first. The same realization problem may arise from 
feedback compensators and even pure derivative feedback. 



3-2 Linear Dependence, 

If the variable parameters are linearly dependent in the coefficients 
of the controlled characteristic equation* then the compensators are 
physically realizable passive or active for every choice of the arbitrary 
roots, because solutions of them are always real numbers. For this reason, 
the variable parameters in the coefficients of the controlled characteris- 
tic equation should be arranged (as possible) such that they are linearly 
dependent. In many cases, the linear dependent relations can be obtained 
by transformations. Consider a single section cascade compensator 



K c (S4 i) 
C5 + PJ 



(3-4) 



If K c , Z and p are variable parameters then the term is a 

quadratic term. However, by defining £ ~0< then (3-4) becomes 

0r c — — ~ttj- f where. K c ==■ °^/ 2 . ^ 3 " 5 ^ 

For the transfer function of (3-5) alone the parameters are linearly de- 
pendent. For the same reason a single section compensator in the feedback 
path the relations can be linearized in the same way. 

The linear dependence of the variable parameter is a sufficient condi- 
tion for the compensator to be realizable* but it is not necessary. The 
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necessary and sufficient condition can be found. For example^ a single 
section cascade compensator the roots within the single section relocation 
zones formed by the rules of reference (7) are the necessary and sufficient 
condition. 

Realization problems may even arise for pure derivative feedback. 
Consider Fig, 3-1 in which the point H X M is assumed to be able to insert 
a summer and the signal at which can be measured. Assumed pure deriva- 
tive feedback is devised as shown in the Fig, 3-1, From rules of signal 
flow graph the determinate of the system contains a term which is the pro- 
duct of the two nontouching loop transmissions. This term contains the 



"X" 




Fig. 3-1. Pure Derivative Feedback 

product of h^ and h^- For a certain choice of the roots, h^ and h^ may be 
complex numbers. For some purposes in control systems, the system may be 
expected to be overdamped and to operate in a given set of roots such as 
discontinuous operation in reference 9. For that purpose the scheme of 
Fig. 3-1 should be avoided. For pure derivative feedback to achieve a 
given set of root certain rules and methods have been worked out. They are 
discussed in detail in Appendix I. 
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3-3 Choice of Arbitrary Roots. 

Any root may be chosen as the arbitrary root. For convenience 
the dominant roots are usually chosen as the arbitrary roots. If there 



dependent variable. If the arbitrary roots are more than two s two of 
them are chosen as the dominant roots , others may be either fixed at a 



roots. Suppose the original system has a zero which is to the right of 
the expected dominant roots. Then in the case of three arbitrary roots 
one of them may be chosen to cancel that zero or be close to it. In 
general, the more the arbitrary roots, the more the design freedom one 
has. The forward gain K may be always treated as a variable parameter, 
so one additional arbitrary root may be introduced even though K is us- 
usally specified within some limits. By treating K as a variable para- 
meter, sometimes it may simplify the relations, but sometimes may make 
the problem more complicated. It depends on the specifications and the 
nature of the compensator, but the increase of design freedom is always 
true. 

For the case when a large number of arbitrary roots i.e., a large 
number of variable parameters is involved, a clear picture of the closed 
loop pole-zero configuration must be assumed. With that pole-zero configura- 
tion the arbitrary roots can be properly related to the dominant roots. In 
the mapping process from the contour of the arbitrary roots to the vari- 
able parameters, only one independent variable is varied, i.e., the para- 
metrical line is obtained in the parametrical space (or coefficient space 
of the reduced characteristic equation). There are infinitely many para- 
metrical lines in the parametrical space, but only those which are within 



is only one arbitrary root 




is the convenient in- 



certain value or be defined to have a certain relation with the dominant 
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the specification limits are of interest* 



3-4 Coefficient as Variables: Mitrovic’s Method* 

In case of coefficients as variables* this method provides a relatively 
simple relation. For lower order systems and for derivative feedback the 
application is direct. For the case of any two coefficients taken as 
variables, the function of the variable coefficients are the same as Mitrovic’s 
coefficient functior which has a certain interpretation from the mapping pro- 
cess. In order to show the coefficient functions are the same as that of 
Mitrovic's method, take a 5th order equation with and as variables 
from Table 1-4, and make substitution of S^+$2 = and ~ n 

from which one obtains: 

O | — Bj- T ~~ 

— Bj - 6 ±u) n + ( - 0 tof 

C Q = 0 Z - - [*{(^-0-*! J 60 * 

fo= ] (3-6) 

f ( = zfoJ n [ 

= + c ^ 4 i s 'V 0 * f wf-Otf (3 _ 7) 

Equations (3-6) and (3-7) are the same expressions as obtained from Mitrovic's 
2 

method . By the same manipulations, functions from Table 1-1 to Table 1-3 
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are the same expressions as obtained from the modified Nitric 0 s method . 
Mitrovic's original coefficient function is derived by mapping a constant 
^ contour through the characteristic equation. By the same mapping 
process* other contours also can be used. In this analytical method the 
functions are expressed in terms of complex variables* the functions for 
any contour can be obtained readily. Because of this extension of Mitrovicfe 
Method to the contour other than constant ^ line* some of the features 
are discussed briefly in the next three sections. 

3-5 Cauchy's Theorem Applied to the Function of Variable Coefficients. 

In the case of two variable coefficients of a characteristic equation* 
the interpretation of the parametrical line of vs f- - is discussed 
in detail in reference (2) and reference (4). Here only the cases of 
variable coefficients which are not in sequence are discussed. Take the 
case of -f-, and f 3 as variable coefficients. Assume the original character- 
istic equation is 

f = + B 3 W B,S + (3-8) 

and assume the controlled characteristic equation is 

f — 5 n 4 - - -- 4 tsS* + B a S*4 f ,S4 B 0 (3-9) 

If -f, and f 3 are so chosen such that f - 0 S then 

F (8*- fj ) S 3 4 ( 6, - f,)s (3-10) 

Assume the arbitrary roots are on a constant J 9 line in the left half of 
the complex plane, then the angle of the S term of (3-10) is / and 
the angle of the S term is /@ if f" 3 and > -f t respectively. 
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If and ‘f, „ then the angle of the corresponding term 

is /-j 0 and respectively as shown in Fig. 3-2(b). The functions 

of -f , and f ^ have negative values as U>« I from Table 1-3. Assume all 

the roots of equation (3-9) are to the left of the f line, then the mapping 
of UJ n from 0 Oo line, must give counterclockwise loops. At the 4 points 
on p(sj curve as shown by (1), (2) (3) and (4), f t and -f must have 
the following values: 

fi < E>i 

-f,= B, (3-11) 

f, ^ B. 

■f, = B, 

By tracing the curve vs -f by (3-11), the vs f 3 curve must have 
a counterclockwise loop and enclose the point ML In general, the curve 
vs f% in which + 1, then a counterclockwise loop of the curve vs / 
implies all the roots are to the left of the ^ line^ i : point M is 
inside the loop. In addition to the conclusions from the reference (2) and 
reference (4), the interpretations of -f , vs f- curve are summarized as 
follows : 



at point (1) 
at point (2) 
at point (3) 
at point (4) 



t,= 

fj = 
f 3 < Bj 



(1) If j = i + 1 ( j = 0„ l s 2,..,), clockwise loop of -f . vs f- 
curve implies all the roots are to the left of the ^ line providing point 
M( is inside the loop. 

(2) If j ^ i+1 ( j = 0, 1, 2,...), counterclockwise loop of 
vs fj curve implies all the roots are to the left of the ^ line pro- 
viding point M is inside the loop. 

The above rules can also be applied when the mapping contour is a 
constant vertical line ( constant). 
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(4 ) 5 - plane 











(c; 



F(s) 

p/ane 



;f , vers f 3 
plane 



Fig. 3-2 Mapping of constant $ contour to coefficient 

plane 
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3-6 Evaluation of roots from the f\ vs ' plot* 

The fx vs fj plot is a mapping process from the complex plane. The 
mapping functions of this process are the functions of the variable co- 
efficients, namely -f. = -j~. (s). In this mapping process regions map 

into region, curves map into curves. An unbounded region on the S-plane 
maps an unbounded region unto the coefficient plane, while bounded regions 
map unto bounded regions. Therefore, the stable root regions discussed in 
the last Chapter map a stable region in the coefficient plane. The map- 
ping process from S-plane to coefficient plane is unique, but the inverse 
(from coefficient plane to S-plane) mapping is not necessarily unique, it 
may have multi-values. By using this multi-value property of inverse map- 
ping, more than one pair of complex roots can be evaluated from the coef- 
ficient plot. Figure 3-3 is a 4th order equation with i~ t and “f ? as vari- 
ables. Fig. 3-3(a) shows the stable root-region and Fig. 3-3(b) shows the 
-f ( vs. f plot. It can be seen from Fig. 3-3 that the unbounded stable 
root region on S-plane map an unbounded region unto f f plane. 

The ^ =0.5 lines on S-plane map a clockwise loop unto the plane. 

Every point on S-plane map a point on plane. However, the inverse 
mapping is double-valued, such as points (2) and (3) on -v f plane 
map two points (2)’ and two points (3)“ on S-plane. Therefore, if suf- 
ficient curves have been drawn on "f\ ^ f plane, by using different con- 
tours (such as constant $ and constant lines), two pairs of complex- 

roots can be evaluated from the f) ^ fj plane directly. For the case of a 
4th order equation, all the roots can be seen on the f f ^ fj plane for 
every choice of f; and f~ • . Real roots are not shown on the coeffici- 

ent plot, but they can be evaluated easily from the reduced characteristic 
equation. 
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Since for every region on the S-plane, there is a corresponding 
region on the coefficient plane, the required sys^-m root-region can be 
mapped on the coefficient plane if such a regv'n exists. For example in 

Fig. 3-3, if the roots are required such roots have 

■// 

and time constant less than 1, then thf root region of the equation has 



a corresponding region on the 






yy 



,~.*y plot which is the intersection re- 

?■ X 



gion of the clockwise loops = 1 and f = 0.5. This region is shown 



yy 



by the shaded area For- > ^ £dier requirements, the region on the coefficient- 
plane which satisfies the requirement can be located by the same procedure. 



// 



Figure 3-4 Figure 3-5 are numerical examples with variables t, 

// 

and t. for 4tV order and 5th order equations respectively. In Figure 3-4 

‘ F 

only one cAista nt Cv_ (= 1) contour is drawn and in Figure 3-5, only one 

/ 

constant* J ( » 0.5) is drawn. From those figures, it can be seen that the 
constant and constant ck. contour constitute curvilinear coordinates 
on the fj f plane. From the point of view of curvilinear coordinates, 
two pairs of complex- roots can be read for every point on the f; f - plane. 



3-7 Stability boundary on the coefficient-plane when the variable coefficients 
are not in sequence: 

As discussed in Section Three of Chapter Two, when the variable 
coefficients are not in sequence, the functions of the variable coefficients 
are undefined on the j <*-> axis. Consequently, the stability of the bound- 
ary line on the coefficient-plane cannot be obtained by mapping the o^axis 
of the S-plane. However, the boundary on the S-plane is well defined as 
has been shown in Chapter Two, therefore, the boundary on the coefficient- 
plane can be obtained by mapping those boundary lines from the S-plane to 
the coefficient plane. This process can be carried out with point by p oint 
calculations or by approximating the boundary on the S-plane by a function. 



80 



The boundary in Fig. 2-12 has been calculated point by point and plotted 
in Fig. 3-6. In Fig. 3-6, a constant ^ =0.1, line is also mapped onto 

the -f ( 'V ^ plane. The ^ = 0.1 line shows a counterclockwise loop which 
implies that there exists a region in which all roots have f >0.1. 
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Fig. 3-4(b) Root-region 
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PART II DESIGN 



CHAPTER IV DESIGN PROCEDURE 

4-1 Specifications . 

The transient response of a control system is determined completely 

if the closed-loop pole-zero (root-zero) configuration is known* The 

contributions to the transient response by each pole and zero can be 

calculated precisely. If a pair of complex roots is dominant , then 

approximations can be made. Based on these dominant roots, the effect of 

all the zeros and all the secondary roots to the transient response can 

6 1 

be visualized from the complete root-zero configuration. s With this 
relation the specifications from the desired response can be conveyed to 
root-zero configuration. As in this method, the design is based on root- 
zero configuration, all specifications must be interpreted in a set of 
appropriate configurations of roots and zeros. Frequency response is 
readily interpreted if dominant roots are assumed. The steady state 
specifications (stiffness and error constant) are contained implicitly in 
the equations and formulas to work with, then they are satisfied automatic 
cally if the final results satisfy all other specifications. However, to 
meet all the specifications may not be possible until some adjustments are 
made. For example, to design a system of high velocity constant (K ) a 
dipole near the origin may be required. If the original dominant roots 
are defined so that all secondary roots are far away to the left of the 
dominant roots, it may result that no root-zero configurations for this 
definition of dominant roots are able to meet the high requirement . 

In this case, the definitions of dominant roots may necessarily be adjust- 
ed so that a dipole to the right of the dominant roots is allowed. 
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In other words, a term of long tail of relatively small magnitude in 
the transient response is allowed. Usually, this long tail term in 
transient response is avoided. But in order to achieve high velocity 
error constant by a relatively simple network, this adjustment may be 
acceptable in consideration of all the factors in design, 

4-2 Design Procedure -- General. 

In the design of control systems, many methods and tools are avail- 
able such as the root-locus method, frequency domain method, pole zero 
configurations, etc. Problems frequently rise about what procedures and 
patterns are more logical. In pure synthesis, the designer is able to 
take the specifications and proceeds logically in a step by step manner 
to a system which will meet these specifications. In these senses many of 
the design techniques are not true synthesis procedures. On the other hand, 
the design procedure proposed by Guillemin* which takes the form analogue 
to modern filter theory is a synthesis technique. In Guillemin’s procedures 
the closed loop transfer function is first determined from t he specifications 
and the compensation is determined subsequently. The design procedure by 
the analytic method developed in this paper is more or less close to Guillemin’s 
procedure. Briefly, the procedures are as follows: 

(1) The regions of closed loop pole-zero configurations are deter- 
mined from the specifications. 

(2) The type of compensator (cascade, feedback^pure derivative, 
or mixed network) is chosen from the consideration of specifications and 
other restrictions. 

(3) Analytic expressions of all the roots and the parameters of 
compensation in terms of the arbitrary roots are obtained, 

(4) Within the pole-zero region defined in (1), choose a set of 
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appropriate pole-zero configurations. Then the parameters of the compen- 
sator are calculated from the expressions in (3). 

Procedure (1) is an interpretation or transition from the specifica- 
tions to the closed loop pole-zero configurations. Procedure (2) depends 
upon the specifications, the nature of the plant and also on engineering 
judgment and experience. These two procedures are general to all methods, 
if the design is based on transient response or time domain. . They are the 
prerequisites. Procedures (3) and (4) are concerned with the method it- 
self. For this analytic method, they are explained in detail in the 
following sections. 

4-3 Design Procedure -- Detail. 

When the type of compensation has been chosen, the design procedure 
for this methid is as follows: 

(1) Formulation of the Characteristic Equation . 

The characteristic equation with the compensator is formu- 
lated in the way such that the coefficients of this equation are functions 
of the free parameters of the compensator. In other words, all parameters 
are independent from one another. 

(2) Formulation of the Root-Coefficient Relations . 

Set the coefficients in (1) equal to the corresponding root 
coefficient. If the order of the equation is "n" there are n n" relations. 
If there are "r" free parameters in the compensator, then ' r" roots are 
arbitrary and the reduced characteristic equation is of the order (nr) 
and has (n-r) coefficients. The root-coefficients relations are formulat- 
ed in such a way that the coefficients are functions of the arbitrary roots 
and the coefficients of the reduced characteristic equation. 
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(3) Formulation of the Functions of the Dependent Variables , 



In the "n" root coefficients relations formulated in Step 2, 
there are "n" dependent variables and r independent variables. The 'W 
dependent variables are the "r" parameters of the compensator and the (n-r) 
coefficients of the reduced characteristic equation. The independent vari- 
ables are the "r" arbitrary chosen roots. Since there are "n" equations and 
"n" dependent variables, the dependent variables can be solved and expressed 
in functions of the arbitrary roots. 

(4) Determination of Stability and Dominance . 

From the reduced characteristic equation, determine the 
stability root-region and the dominant root-region. The dominant root 
region depends upon the degree and the definition of the dominance which 
is defined by the designer from the consideration of specifications. 

(5) Evaluation of the Compensator Parameters. 

If the desired dominant roots are within the dominant root- 
region which has been determined in step (4), then choose the dominant 
roots (usually independent variables) and substitute them into the formula 
of step (3), then the parameters of the compensator are evaluated. 

(6) Realization of the Compensator Network . 

If the parameters of the compensator evaluated in the last 
step are within the specified range and can be realized, then the design has 
been completed. 

(7) Design Adjustment . 

If the value of the compensator parameters are not in the 
specified range or the desired dominant roots are not in the dominant-root 
region, then it indicates either this type of compensation cannot meet the 
requirements or the definition of the dominant roots must be adjusted. When 
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other types of compensation have been chosen or the degree of dominance 
has been adjusted, the design procedure is just a repeat of the above 
steps. 

The procedures described above are general. The first three steps 
are the partitioning process as discussed in Chapter I. For some specific 
cases, formula have been derived as shown in the subsequent sections, then 
the design is just a computation and is started from Step (4). The compu- 
tation is best carried out by tabulating all previous data and by pro- 
gressing computation. Once those data are tabulated, all information 
can be found from the table and if some changes of the type of the com- 
pensation are made, most of the data can be utilized. 

4-4 Free Parameters of the Compensator. 

In step (2) of the last section, a terminology ,f Free Parameters" was 
used. It may need further explanation. Consider a cascade single section 

R-C compensator. Let the transfer function be 

K r (S + Z) 

S + P 

The apparent parameters are K c , Z and j? . Assume the system is type I 
and the velocity error constant is specified. If the forward gain K 
of the uncompensated system was determined by the requirement and in 
the partitioning processes, K was treated as a given numerical constant 
then K c is not arbitrary in order to meet the requirement of K^. In other 
words K c must satisfy the following relation: 

K c = P/z 

Thus the three apparent compensator parameters K c , p , and Z are related 
by the requirement of K^. Only two of them are independent variables, the 
third one is dependent. The independent variables of the compensator in 
this sense are given the terminology "free parameters". 
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For other cases, the same treatment can be applied to the compen- 
sator parameters. In the above example, if the forward gain K was also 
treated as an adjustable parameter, then all the three apparent para- 
meters are independent because K c K is also independent. However, the re- 
sult of the partitioning processes are different for the two cases. But 
the design can be carried out by either case. In the former case (K is 
constant) there are two arbitrary roots, while in the latter case (K is ad- 
justable) there are three arbitrary roots. The expressions of the para- 
meters are different for different cases, but the design can be carried 
out in both cases. 
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CHAPTER V DESIGN OF PURE DERIVATIVE FEEDBACK 



5-1 Tachometer Feedback - Second Order. 

For compensating lightly damped instrument servos, tachometer feed- 
back is commonly used. For second order cases, as all information can 
be determined precisely, only the formulas are presented. 




Fig. 5-1 Second Order System with Tachometer 

Feedback 

The block diagram of a Type 1 second order system is shown in Fig. 
5-1. The separation of the forward gain by and is for the purpose 
of proper adjustment of h. Assume the specifications are given and 

^ . The design procedures are as follows: 

(1) Formulation of the Characteristic Equation 

The Characteristic Equation is 

5^+ ( 0.4 K*h )S 4 K,K 2 = 0 (5-1) 

(2) Formulation of the Root-Coefficient Relation. 

Since there is only one parameter of compensation, there is only one 
independent variable. The independent variable is not necessarily the 
arbitrary root itself, it can be a function of the roots. In the case 
of one arbitrary root and that root has to be complex, the independent 
variable can be chosen as either f or CO h of the pair of complex 
conjugate roots. For second order case, the only choice is ? , 
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because CO n is fixed by the gain. By doing this the root-coefficient 
relations are: 

a -I- K a h — t (5=2-1) 

K,K Z = ( 5 - 2 - 2 ) 

(3) Formulation of the Functions of the Dependent Variables 0 
Since the root of the reduced characteristic equation is the domin- 
ant root itself, there are no stability and dominant region problems. 

Here the dependent variables are cG n and h( ^ is chosen as independent 
variable). Their expressions as functions of the independent variable are 
as follows: 

From (5-2-2): 

w n = JZkI (5-3-1) 

From (5-2-1): K^h = 2JW„ - <X = If - CL 

iS/fTST-O- (5-3-2) 

H = — 

For the given value of J , h is computed. 

5-2 Third Order: Tachometer Feedback. 

The block diagram of a third order type one system with tachometer 
feedback all around the plant is shown in Fig. 5-2. 



"0 1 cy 



5(5taJCS tb) 



hS 



Fig. 5-2 Third Order System with Tachometer 
Feedback 
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(5-3) 



The characteristic equation of the uncompensated system is 

5 3 + B 2 s 2 + b,s -f B 0 =o 

where B 2 = a 4* b, B^ = ab, B q = 

Assume the forward gain K and the damping ratio f based on the 
dominant roots are given. The design procedures are: 

(1) Formulation of the characteristic equation: It is easily 

formulated as follows: 

5*+ B 2 S 2 + ( B , + KJiJs -f B 0 — O < 5 * 4 > 

(2) Formulation of the root-coefficient relation: Since 

there is only one control parameter h, only one root is arbitrary. 

This root can be chosen either the real root or ^ of the complex roots. 
Assume the real root is chosen and designated as (3^ is positive on 

the negative real axis). Then is the independent variable. The re- 
duced characteristic equation is of order 2 and therefore has the form 

5 * + C,S + C a = O 

Since the roots of this equation are actually the dominant roots, it is 
expressed by j* and CxJ n for convenience. Here and U) n are not in- 
dependent variables but dependent variables. By the choice of the 
above definition, the root-coefficients relations are: 

Bi = ti, (5-5-1) 

B, + K d h = (5-5-2) 

B 0 = <^n S . < 5 " 5 “ 3 > 

(3) Determination of the functions of the dependent variables. 

Here the three dependent variables are ^ , u> n and h. In the 

actual computation, only one dependent variable must be found explicitly 
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as a function of the independent variable, the others can be computed 
by using the implicitly functions of formulas (5-5). But for the pur- 
pose of illustration, they are solved and shown in equations (5-6) 



s = (B.-s.jjs '/zjrs: 


(5-6-1) 


— J 


(5-6-2) 


K 2 h = B 0 /s, + B,s, -5,'- B, 


(5-6-3) 



Note the three variables £ , tJ n and h are expressed as functions of the 
independent variable. 

(4) Determination of stability and dominant root region: 

In this one control parameter case, the stability and dominant root 
regions are lines. These lines are the same as the root locus, by treat- 
ing h as variable and K as a fixed number. The only difference between 
them is the independent variable. In the root- locus method, the indepen- 
dent variable is h while in this method it is the arbitrary root S^. A 
numerical example is shown in Fig. 5-4. The degree of dominance can be 
visualized easily from those plots. 

For this particular scheme, it indicates a max. damping. This max. 
damping can be found from equation (5-6-1) analytically because ^ is ex- 
pressed explicitly as a function of S^. Differentiate ^ with respect to 
Sj and set it to zero, one obtains: 

S,„, = b,/j 

Substitute into equation (5-6-1), obtain 

sL = 

Substitute this value of into ^ and C0 n , the real part of the 
complex root becomes 



(5-7) 

(5-8) 
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fu) n — = S, m 



(5*9) 



Equation (A- 7) shows that if the tachometer loop gain I^h is so 

g 

adjusted such that the real root is ( — — - ) then the complex 

roots have max. damping for the given value of B ( * K) and this max. 

o 

damping is shown in equation (5-8). Equation (5-9) shows that for max. 
damping the real part of the complex roots is equal to the real root. 
Fig. 5-3 shows this root configuration. 




Fig. 5-3 Root Configuration for max. f of 3rd Order 
System with Tachometer Feedback. 

For this root configuration, the tachometer loop gain and oJ n are 
obtained from equations (5-6-2) and (5-6-3) as follows: 



iy. = nr 4 K-s, 



(5-10) 



L. ) = 3 B ° 

Wn J 

Equation (5-8) gives a quick look of the max. possible damping for given 
K or the max. K for a given ^ in a third order system with tachometer 



(5-11) 
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Fig. 5-4(a) Plot of Equations (5-6-1) and (5-6-2) for numerical values 

given in example' 5-1. 
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feedback. 



Example 5-1: For the block diagram shown in Fig. 5-2, if b « 10, 

a = 5, = 750; the uncompensated system has a pair of complex roots 

on j 60 axis. For this pole-zero configuration: B 2 = 15, = 50, = 750. 

From equation (5-8), the max. possible damping 



5 = 



_ IS i 



Z?X 7 <rO 



o.q-o '3 



From (5-7) S^ m = 5. 



Assume this value of f is taken, then from (5-6-2) and (5-6-3) 

<^> n = 12.3 



K 2 h = 18.5 

For other root configurations, ^ and OJ n as a function of are plotted 
in Fig. (5-4) (a) and the root-locus for h as variable is shown in Fig. (5-4) 
(b). If the dominant root is defined so that the secondary root has greater 
real part than that of the dominant roots, then the portion of the root 
locus when Sj 5 5 is the dominant root region. 



5-3 Third Order: Tachometer Feedback to Interior Node. 

Consider Fig. 5-5 in which point "X" is assumed to be available to 
insert a summer. The analytic design procedures are; 




Fig. 5-5 Tachometer Feedback: Minor Loop 
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(1) The characteristic equation is: 



5 3 + ( 6 2 + ( B,+ bK z h J-S + & o —0 (5-12) 

where B~ = a+b, B. » ab, B = K,K- 
2 ’1 o 1 2 

(2) Root-coefficient relationship. For the same definition 
of dependent and independent variables as described in the last section, 
the relations are 



8 x i K z h = *?u) n + 5, 

B, + bk,h = + Zju; n S, 

B 0 — 5 i 

(3) Function of dependent variables ( 

From equation (5-13-1) and (5-13-2), eliminate K^h 

J== ( ?60 n (b-S,J 

Substitute C*J n from (5-13-3) into the above expression obtain 



OJ 



(5-13-1) 

(5-13-2) 

(5-13-3) 



K*h ). 






B,b -6, + -|? - bS, 



2 ( b-Sj J 8b /s, 

The other two variables can be obtained in the same way but are not 



(5-14) 



necessary. 



(4) Stability and dominant root region: For the numerical 
values given in Example 5-1, the dominant roots ^ and cO^ are plotted 
in Fig. 5-6 as functions of from equation (5-14). The degree of dom- 
inance can be seen easily. 

(5) Take = 0.5 which corresponding to = 20. From (5-13- 



3) 



0 CO n = 6.14 

From either (5-13-1) or (5-13-2), obtain 



Kzh = 11.14 
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_ 5 , ( real rooi) 



o 



T 



T 



o 



\_^o O 



Fig. 5-6 f and (jJ n as a function of the real 
roDt of the systems shown in Fig. 5-5. 




100 



5-4 First and Second Derivative Feedback, 



In case of first and second derivative feedback around the plant 
for a given value of forward gain, only the coefficients of the first 
and second derivative term are variable. The design may be carried out 
using two variables. But for more freedom it is better to take the con- 
stant term also as a variable then three arbitrary roots are available. 

The coefficients of the reduced characteristic equation and functions of 
the variable coefficients have been derived and tabulated in Table 1-6, 
Example 5-2. Given the original characteristic equation, 

S>* + iOS i 4 3IS 2 4-303 + K ~ O (5-15) 

The forward gain K must be greater than 1000. The gain at the stability 

limit for the uncompensated system is 84, first and second derivative feed- 
back are chosen. When three variable coefficients are used, equation (5- 
15) becomes 

5% |03 3 4 f 2 3* -f f , 5 + f 0 ~ O (5-16) 

where 4- = K 

1 o 

From Table 1-6, the coefficient of the reduced characteristic equa- 



tion and the functions of the variable coefficients are: 

c c = e^-c-s.+^ + Sij ( 5 - 17 -D 

f 0 = CoS.SiSj (5-17-2) 

f , = c„ ( 5,5, -1 ^ S, J (5-17-3) 

f a = C 0 (S,-t5, + S 3 j +(J A +i 2 S J +J J 5,j (5-17-4) 



where (-S^), (- S 2 ) and (-S^) are the arbitrary roots and C q = S^. Here 
the order of the equation is four, all the roots may be complex. There 
is only one constrained root so it is convenient to replace by it- 
self. The order of the reduced characteristic equation is one; equation 
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(5-17-1) alone determines all the roots. Assume (-S^) and (-S 2 ) are the 
dominant roots, then (5-17-1) becomes 



10 = + S, +i + (5-18-1) 

f„=^S J S 4 (5-18-2) 

There are no zeros in this system, all secondary roots (-S^ and -S^) are 
required to be to the left of the dominant roots. Assume (-S^) and (-S^) 
are real roots, then for the definition of dominance, the max. K occurs 
when . From (5-18-1) for this value of the secondary 



roots, the conditions are: 



fu) )y = Z.5 
max. K = f 






(5-19-1) 

(5-19-2) 



For ^ =0.5, the max. K = 156. This is far below the required K, then 

/ , 

the secondary roots must be complex. Assume J and uJ n denote the second 
roots, then equations (5-18) become 



(5-20-1) 

(5-20-2) 



4 ftd n ~ 5 

K = tol uJ? 

Notice in equation (5-20-1) there are three arbitrary variables and this 
equation alone determines all the roots. Assume the dominance is defined 
as ^ tO n then from (5-20-1), the dominant root region is 

< 2.5, which is shown in Fig. 5-7a. Within this region choose 
^ =0.5, = 2, then from (5-20-1) the condition for secondary 

t / 

roots is ^ (-0,, = 4 

Substitute into (5-20-2), obtain j<^ — (o ^\- / ^ > ' i 
For K = 1000, obtain ^ ~ 0.253. For this choice of the dominant roots 
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and f 2 are readily calculated from (5-17-3) and (5-17-4) as follows 

f. = 532 
f 2 = 270 

The root configuration is shown in Fig. 5-7b. For another choice of the 
dominant roots, the computation can be carried out in just the same way. 
Example 5-3: Given the original characteristic equation 

5 r 4 (28^4 £-205*4 l5oo. S 4 K —O 

The forward gain K must be greater than 2000. For this value of K the 
system is badly damped, first and second derivative feedback are chosen. 
By using 




Fig. 5-7a 

Dominant Root Region of 
Example 5-2 



Fig. 5-7b 

Root -configuration of Example 5-2 
for = 0.5, (J n * 2 
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three variable coefficients, the working formulas from Table 1-6 are as 
follows : 



C,= 


S h -CS, + 5 t 4S 3 ) 


(5-2i- n 


Co = 


B 3 -C, (5, + S i +S 3 j_(^+^S^ + 5 J S ( ) 


(5-21-2) 


II 


l,S*SsCe 


(5-21-3) 


f . = 


C 0 ( -5i S 2 4 4 .Sj S ( J 4 C, 


(5-21-4) 


II 


Co + C, ( S,5 Z + 5^ t S,5,j) 4 S,5,5 J 


(5-21-5) 



Assume (-S^) and ( - S 2 ) are the dominant roots. The order of the equation is 
five, one root must be real. Let this real root be (-S^). Here three roots 
are arbitrary. may assume any value. For dominance, assume 



5^ = ajo On (5-22) 

Substitute (5-22), B^ = 16 and = 128 into equation (5-21), one obtains: 

C,= IC? - 4 (5-23-1) 

Co= US -( l + (5-23-2) 

f„ = t?Lo:,c„ < 5 " 23 - 3 > 

f, = C.lJ*(l+«5 (5-23-4) 

= 4fU„C. <-(1+ Sy^C, +2fUJ 3 (5-23-5) 



Notice (5-23-1) and (5-23-2) determine all the roots, (5-23-3) determines 
the forward given K while (5-23-4) and (5-23-5) determine the derivative 
gain. Since there are no zeros in the system, the secondary roots are 
required to be to the feft of the dominant roots. Let a = CO n , the 
transformed coefficients from Table 2-4 are as follows: 
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(5-24-1) 



D, — C,~ 

T> 0 = C 0 (5-24-2) 

The stability boundary lines are determined by C^^and C Q ^>0 . 

The dominant boundary lines are determined by J) >o and J)^ > O 

They are shown in Fig. 5-8a. For ^ =0.5, the results of computation of 

C , C-, f , D and D- are tabulated in Table 5-1. Examine Table 5-1, 
o* V o’ o 1 * 

when 0O n is between 4 and 5 the gain K is greater than the specified 
value 2000. This range of (jJ n is within the dominant region since both 
and D q are positive. The root-configuration for U) n = 4, f = 0.5, 
is shown in Fig. 5-8b. For this choice of the dominant roots, the vari- 
able coefficients are computed readily from (5-23-3), (5-23-4) and (5-23-5). 

f = K = 2040 
o 

fj = 1534 (original 1300) 

= 576 (original 520) 

For other values of and other choices of the dominant root the 
computation is the same. 



UJ n 


c, 


Co 


fo 




Po 


2 


12 


72 


576 


+ 


+ 


3 


10 


50 


1350 


+ 


+ 


4 


8 


32 


2040 


+ 


+ 


5 


6 


18 


2250 


+ 


+ 


6 


4 


8 


1728 


- 


+ 


7 


2 


2 


686 




+ 


8 


0 


0 


0 







Table 5-1 Example 5-3 for ^ = 0.5 
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boundary of 
stability 



\ 



-CK. '4 



fW,=: 2.66 

I 

k 
k 
k 
tf 

I 

fc 



boundary 

of 

dominance 

Ol- 



-3 "2 



J CO 



3 

z 

I 



Fig. 5-8a Stability and Dominant Root Regions of Example 5-3 



l 



Fig. 




5-8b 



Root Configuration of Example 5-3 for 
U) n = 4. 






0.5 
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CHAPTER VI 



FEEDBACK COMPENSATOR DESIGN 

6-1 Introduction. 

The compensation by pure derivative feedback discussed in the last 
chapter usually adjusts only some of the coefficients of the character- 
istic equation. The effectiveness of this type of compensation is less 
than that if the compensator parameters enter all the coefficients. This 
has been shown by the reduced characteristic equation in Chapter II. 
Moreover, the tachometer feedback reduced K and acceleration feedback re- 
duced and so on; and a system can seldom be stabilized without adjust- 
ing the lower order coefficients of the characteristic equation. There- 
fore, if the compensator is chosen in the feedback path, the compensator 
often consists of a tachometer followed by a network. The tachometer 
effectively changes the coefficient of the first derivative term, while 
the poles and zeros of the network enter all the coefficients in the way 
of the coherent nature of algebraic equations. This chapter is a discus- 
sion about the design of such a compensator by the analytic method. 

6-2 Second Order Compensator. 

Consider Fig. 6-1 which is a second order system with feedback compen- 
sator consisting of a tachometer followed by a lead network. 




Fig. 6-1 Second Order System with Feedback Compensator. 
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The characteristic equation of the compensated system is 

_5%( B. + p + ^4*. ) S + ( 6 C + B, pJ ^ ■+ &ol° — O (6-i) 

where = a, B q = In this characteristic equation* there are 

two variables in the coefficients, therefore two roots are arbitrary 

if and p have no restrictions. The reduced characteristic equation 

is of order one and has the form of S + C =0. The root coefficient re- 

o 

lations are: 



B, + p + K z = 2^iO n + Co 


(6-2-1) 


Bo + &,jp — <*>n + 


(6-2-2) 


6> 0 tp 


*,6-2-3) 



In equation (6-2), the dependent variables are jp , and C^; the in- 
dependent variables are ^ and 60^ are assumed to define the 

dominant roots. determines the stability and dominant root region* 

while p and must assume values according to the dominant roots. 

Since equations (6-2) are linear functions of the dependent variables* 
it is convenient to write in the matrix form as follows: 



/ 


1 


- 1 








B,“ 


0 


B, 






P 




^ - B„ 


1 


B e 






1 

o 

o 

1 




0 



The three dependent variables can be solved readily, but only 
expressed as an explicit function of the independent variables, 
rule, one obtains ? 



c 



o 






must be 
By Cramer's 



(6-4) 
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C in this case is the actual real root. The stability and dominant 
o 

root region are determined by equation (6-4) by varying ^ and „ 

For design purposes, only those values of ^ and which are within 

the specifications are computed. When ^ and have been chosen, 

^ka and p are calculated and the design is completed! 

Example 6-1: In the block diagram Fig. 6-1, = » a = 23. 

Therefore B = 10“* and Bj = 23. From equation (6-4): 

0 ax(o 5 J u> n - 



The stability criterion is 



I0 S - ^Jn 



zxio r % uJ n -z5uJ* 



O l.e, 

>0 



or 



10 s - cJ,* > <0 and 



ZXIO fOJ n - o 

These two inequalities define the stability region and is shown in Fig. 
6-2. If the dominant region is defined by C 0 ^>2.f^J n then the 



criterion for this definition is 

IO S - 



JI- > 



2X10*?^- 

This region is also shown in Fig. 6-2. For ^ =0.5, the range of C d n is 
0<C U) n -C 220. The values of C c , p and ^ka for the dominant roots in 
this region ( ^ = 0.5) are computed as shown in Table 6-1. The final 

decision of the design is to pick up a value which is favorable to all 



compensator parameters and other specifications. From Table 6-1 assume: 





Co 


P 


Ki&a. 


10 


io 4 


10 


9770 


20 


5000 


20 


4960 


40 


2480 


38.4 


2460 


60 


1600 


56 


1580 


80 


1200 


77 


1180 


100 


922 


92 


900 


200 


312 


125 


365 



Table 6-1 Values of jo , l^ka and C q 

of Example 5-1 for J =0.5 
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Fig. 6-2 Stable root-region and Dominant root-region 

of Example 6-1 



i 



110 



The dominant roots are chosen to be ^ = 0.5, Qj = 100 s then the third 
real root is -922, p = 92 and = 900. The closed loop transfer func- 

tion is 



0c 

0 o 



K, KJ^ P) 



(S + S,)fS + 50(s + 53 ) 



The closed loop pole-zero configuration of the above choice of the domin- 
ant roots is shown in Fig, 6-3* 




Fig, 6-3 Closed Loop Pole-zero Configuration of 
Example 6-1 for the dominant roots 
J = 0.5, U) n = 100 

From Fig. 6-3, it can be seen that the system is essentially second 
order because the real root is far away from the origin. The effect of 
p to the transient response is also easily visualized from the pole zero 
configuration. 

6-3 Third Order. Compensator. 

Fig. 6-4 is a block diagram of a third order system with feedback 
compensator. Assume the forward gain is fixed from the steady state 

accuracy and the velocity constant must be greater than a given value K . 



Ill 




Fig. 6-4 Third Order System with Feedback Compensator. 
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The characteristic equation of the compensated system is 



5 f + ( S £ + PJ5 3 + ( B, + Bp * KjJ^+( S,p* Kj t neJS i d„p=0 (6-5) 



where — &- + ^ J 3 ~ & b , B 0 ~ K.K* 

In the above equation, there are three variables, namely fp ^ , 
However, since is specified, those three variables are not 
one to another, they are constrained by the specification 
compensated system is 



jo and £ 
independent 

K of the 

v 



K 



V 






(6-6) 



Rearrange equation (6-6) 

K 2 M -PB, 

Substitute this expression into the coefficient of the first derivative 
term of equation (6-5), it becomes 

B,p + K £ t l t B 0 = B,f> + -pB, ♦ B, = £jh + Bo 

Then equation (6-5) becomes 



S 4 +(B i +pJS i +(B,tB 1 p + lf^ t )5 2 +(-^P + B 0 )S 4 B o p = 0 < 6 ' 7 > 

The coefficients of characteristic equation now have only two variables*, 
namely £ and jD . The specification of eliminates one of the vari- 
ables of the compensator. This is always true. The more specifications 
given, the less the freedom of the compensator. Since there are two vari- 
ables in the coefficients, two roots are arbitrary. The reduced character- 
istic equation is of order two and has the form: 
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B, +t> = + C, 



( 6 - 8 - 1 ) 



6,+ B,p+ = UJ* 4 C 0 4 2fW w C, 

Be 



jp + B 0 = UJ 0 C, 4 *$U) n C t 



K„ 

E>P = <C 0 



( 6 - 8 - 2 ) 

(6-8-3) 

(6-8-4) 



where ^ and 60 n are the arbitrary roots while C q and are the co- 
efficient of the reduced characteristic equation. As the four dependent 
variables C-, C , p and in the above equations are linear functions 8 

it is convenient to express by matrix form as follows: 



O 



O 



~w! 



B. 



-2ju) n 

-u) z 



0 

0 





r C," 








Co 




^ - B, 




1 0 




B 0 




kJt 




0 



(6-9) 



By Cramer's rule, one obtains 

C = 



B*- 

MW- Y"n * ^ 



(6-10) 



The stability and dominant region are determined by C q and C^. Choose 
dominant roots the parameters > p and £. of the compensator are 

calculated from equations (6-8) and (6-6) 

Example 6-2. For the system of Fig. 6-4 given a = 10, b = 30, = 

21,000^ K v =r 60, damping ratio ^ of the dominant roots be approximately 
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0.7 and C O n = 10. 



For the given values, = 40 



Bj = 300 

B = 21,000 
o 



Substitute into equation (6-10) for ^ = 0.7 



S' ,2 



c = 



44 XIO - 8. 4*70 <^n 4 1-4 u) n 
ajf* (0 4 td„ - 35XlO l uJ* 4 tJ, 4 



(6-11) 



For (J n = 10, C o = 1321 

From (6-8-4) p = 6.26 

From (6-8-1) C, = \ + P ~ Z$U) n = 32.3 

0^ and C ^ determine the secondary roots, they must be checked before 

going further. The reduced characteristic equation is 

3* + 3Z2 S + 1520 = O 

The roots of this equation are -16.1 + j 32.5. Compare the real part of 
these roots to that of the dominant roots ( = 7)* it indicates the 

presumed dominant roots are dominant. 

From (6-8-2) 

Mt = + c 0 ♦ 2 fu)„c,-B,-e, 2 p — on 

and finally from (6-6) 

-f- = izi, , Z = (4. (5 

The compensator becomes 

I3a3 x CS f 1 4,1$-) _ 

K z ( 5 t <o.z<o) 

The calculations above are for one point only. The choice of the 
dominant roots does not insure a satisfactory result. For more freedom 
and adjustment, the stability and dominant root region are calculated 
first as shown in Fig. 6-5. From Fig. 6-5, adjustments of the dominant 
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Fig. 6- 



Fig. 6-1 

) 
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roots can be made. 



The closed loop transfer function is readily formulated as 
0c _ K,fG(^P ) 

e R — cs+s ( ) cs^-s a )cs + s 3 )(s+s*) 

The closed loop pole-zero configuration for ? = 0.7, uj s 10 

/ n 

(dominant roots) is shown in Fig. (6-6). 

Example: 6-3. Given the same plant and the compensator as example 
6-2 except the specifications of is relaxed. Here are three free ad- 
justable parameters, three arbitrary roots are available. The control- 
led characteristic equation is the same as (6-5), the root-coefficient 
relation are as follows: 



6,t p = 2?aJ n +^+5 4 (6-12-1) 

6,-* %P + ~ + V* (6-12-2) 

B,p + K t g t i t K = + 2f (6-12-3) 

Kp=6j"s3\ (6-12-fj 



The four dependent variables are , Jo , and the three independent 

variables are ^ , U and S^. Assume and are complex, and are 

denoted by f and CJ n y . If (jj' is chosen as the dependent variable, then 
from (6-12-1) and (6-12-4) one obtains 





BJ — O 



(6=13) 



This equation alone determines all the roots since it contains , Cd n , 
, and (jJ^ . Assume ^ (dominant root) =0.5, ^ ( secondary root) 

= 0.7, and substitute the numerical values of K = 21,000, “ 300 and 

B 2 = 40 into (6-13) and solve for Li) n > one obtains 




1^1X10 



4- 

■[>*r 



I + V+t/o-fyCcJ^+o) 



(6=14) 
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Assume only 10 — GiJ ^ ^ 20 are of Interest, the results of compu 

tations from Equations (6-12) and (6-14) are tabulated in Table 6-2. 

Notice in Table 6-3 when UJ n = 10 and UJ n = 23. 5^ Z is negative. 



(/0 n 


Wn 


P 


Kj t 


Z 


?/z 


10 


270 


348 


62660 


10.25 


34 


23.5 


2.63 


575 


-7.25 


-3.62 


15 


109.2 

21.4 


128.5 


9085 


17.2 


7.48 


20 


54.4 


56.4 


5324 


22.3 


2.53 


19.6 


Table 6-2 


Example 6-3 for 


f = 0.5, 


% '= 0.7 




For 


simple 


realization of 


the compensator. 


the roots 


of UJ n = 15 



(jj ' = 109.2 are chosen. For this choice of the roots, the compensator 



is 




(5+ '28.5.) 



where 




Kj_it 



the system are P =0.5, 
is a zero at -128.5. From 
roots ^ =0.5, = 



If Kj = 1, = 0.433. -f ie four roots of 

<jJ = 15, f / = 0.7, id' = 109.2. There 

the closed loop pole-zero configuration, the 
15 have high degree of dominance. 
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CHAPTER VII 



CASCADE COMPENSATION DESIGN 

7-1 Single Section: 

As cascade compensation is commonly used, the design procedure dis- 
cussed in Chapter IV, is applied. In Chapter I a third order system was 
analyzed, here a nth order system in the general case is considered. 



p— 


0((5+i) 




K 








r n i a 

$ -*■ V. s + ~ 




i 

1 












1 





Fig. 7-1 Single Section Cascade Compensation for a 

nth Order System 

Consider the system of Fig. 7-1. Assume the characteristic equa- 
tion of the uncompensated system is 

s n + 6 n _ ( 5 n " < + B n . 2 5 n '% --- + B z S Z + 6,S -h B o =^0 (7-1) 



For specified forward gain K and error constant (K^, K^, etc.) the 
three compensator parameters Z, jp and 0( are not independent one 
to another, but related, namely 




The characteristic equation of the compensated system is 



(7-2) 



If! ' 

-S +(8„+pJS ^ + (8;.+ (7-3) 

+ B „jp=0 
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In equation (7-3), there are two variables jp and o( in the coef- 



ficients, then two roots are arbitrary. The reduced characteristic 
equation is of order (n-1), and has the following form. 



n-i 



S" ' + C„. 2 s" _1 + --- -+ c,s-tC c =o 

If the two arbitrary roots are presumed to be the dominant roots of the 
compensated system and are expressed by j and 
coefficient relations of equation (7-3) are: 

z 



CU ] then the root- 



B<* + 



B, 

B, 



+ 

+ 



B.P 
B iP 

B,p 



= C,w„ 
= 



Q tfw. 



C 2^n + C,ljH 4- Co 

4 C, 



(7-4) 



B. 

A, 



= c 4« + Op* + o 2 



The N+l dependent variables ( p , o( , C q , 2 ) are linear 

functions, then equation (7-4) can be expressed in matrix form 



3 n 0 O O 

Bn-. 0 O 

B„-, O ' 

I I 

B„-j 0 • , 

* \ 1 

, . 

, \ | | 

* I | » 

■ . 1 

, 0 4> 

I 

; 0 o - 1 

B 2 0 - | 

B| 

Bo 0 -u 2 z O 







_ 




1 


0 -1 




P 






- 1 -2fH> 




c< 




< 


”2(K -w* 




C 0 




0 


-wj 1 




c. 




0 


0 * 

I 

1 

1 

1 1 

1 t 




G 

\ 

( 

1 

l 


1 


1 

l 

I 

1 

1 


1 1 

1 1 

1 ' 

i I 

i 1 

O 




1 

1 

* 

_ C "G 


j 

i 

! 

| 

[ 


1 

1 

, 

0 






(7-5) 



I 

♦ 

i 



6, 

O 

0 
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By Cramer's rule one obtains 



-eJ'fiB, + + y,B,u)* 4 - - - 

e„ -t <f, + % e, u)Z + % e,^: + 



-• %B- cj n ~ 



e„M 



/ n 



, + • 



(7-6) 



where c f ? , — — I 

fz - 

,<& = -4f + 8f 

! (7-7) 

l 

I 

I 

The other dependent variables can be expressed in the same way but this 
is not necessary as far as computational purposes are concerned. When C q 
has been calculated, the other variables are readily evaluated from equa- 
tions (7-4). Choose ^ , evaluate the stability and dominant root 

interval of . Within the dominant root region, choose u) n and evalu 

ate p and 0( . If the network is realizable, the design wrok is 

completed. 

Example 7-1. Given Q* = r and K is not to be reduced. Then 
the characteristic equation of the uncompensated system is 

S 2 + 400 = 0 

B- = 1, B = 0, B = 400, B = 0 where i 2 
2 1 o i 

The order of the compensated system is three and that of the reduced 

characteristic equation is one. Then 

C, = 1, C . = 0 when i ^ 1 

1 ’ i 
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From Equation 7-4: 



B 0 p — C c 

4 B,p = C,CJ* 4 C 0 zfU) h 

6, 4 B z p = C^JOOn 4 Co 

Substitute the values of B's into the above equations, obtains 

q. 00 p = Co^n 

^■00 <X. — - 4 ClJpC-o 

p = 2^UJ 0 4 Co 
From equation (7-6) and (7-7) 

c _ 6 .( 

° 6 „- w„‘ 

Substitute the values of B's, obtain 



C 0 = 



Soo fu) n 



GJ* - 400 

For stability the criterion is C Q l> 0> namely 



gOO f OJn 



> O 



CO,, 1 - 400 

The numerator never becomes negative, then the inequality becomes 



U0 L 



4 - 0 O 



Define the dominant roots such that the real part of the secondary 
must be at least twice that of the dominant roots, then 

X oo 



Co ^ / x 

' co* - <U>o 

Then the criterion for dominant roots is 

(loo ^ uJ* c Sov 

Those regions are shown in Fig. 7-2. 



Z f CO?) 



(7-8) 



(7-9) 



(7-10) 



(7-11) 



(7-12) 

roots 
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Within the dominant root region, for j 3 =0.5, and ^ = 0 . 6 ^ the 

compensator parameters are computed from equation (7-8) for different 
values of cO n and tabulated in Table 7-1. Table 7-1 includes not only 
the compensator parameters, but also all the roots of the system for dif- 
ferent choice of the compensator. Assume the dominant roots are chosen 



as 


f = 0.6, 
= 4.41, 


U) =30, then 


from Table 7-1, C = 28.8, 
o 

14.8, The compensator is 


p = 64.8, 
tc ~ 




f — O.S 


f = o, (o 




Co 


P 




Co 


P 




22 


220 


242 


13.2 


264 


291 


14 


24 


54.6 


78.5 


4.72 


67.7 


94.3 


5.34 


26 


37.7 


63.4 


4.12 


45.3 


76.1 


4.6 


28 


29.2 


57.2 


4 


35.1 


68.5 


4.41 


30 


24 


54 


4.05 


28.8 


64.8 


4.41 


32 


20.5 


54.6 


4.2 


24.6 


65.6 


4.53 


34 


19.3 


55.7 


4.53 


23.2 


66.8 


4.87 



Table 7-1 Computation Result of Example 7-1 
The closed loop-pole-zero configuration for this choice of the dominant 
roots is shown in Fig. 7-3(a). In Fig, 7-3(a) it can be visualized that 
the peaking time is small but the overshoot may be too large because the 

g 

zero is too close to the origin . If ^ = 0.6, = 24, are chosen as 

the dominant roots, then from Table 7-1, C q = 67.7, = 94.3 A Of = 5.34 

and £ =17.7. The compensator is 

__ s.wc mill - 
— (s -+ 

The pole-zero configuration for this choice of the dominant roots is shown 
in Fig. 7-3 (b). For this choice of the dominant roots the peaking time is 
still small, but the overshoot is reduced. 
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Fig. 



.. Fig. 



Fig. 




7-2 

Stable and Dominant Region of Example 7-1 




7—5 ( a) 

Closed Loop Pole-zero Configuration for Dominar* 

Roots ^ = 0.6, u) n = 30 for Example 7-1. 




-3(b) 

Closed Loop Pole-zero Configuration for Dominant 
Roots ^ = 0.6, 24 of Example 7-1. 
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^ 300 

Example 7-2; Given Q- = 5(<,+2.)($-t-5j an d *" S n0t to 
reduced. The characteristic equation of the uncompensated system is 

S 3 + 1$ + 10^ + 300 ; 0 

From Equation 7-4: 

7 + p = c, + 

10 4 7P = u)* + C a + 2fU>nC, (7-12) 

lOp + 3000C = C,^ 4 ; 3^0^? 

From Equation (7-6) and (7-7) 



- .1 
C o c0 n 



c. = 



For f =0.5 



For 



= 0.3 



-3oo [ - 1 0 + f4- f U-),i 4 ( i- 4-f ) u) n J 
300 - 7 4 

( 7^ n ~io) 3oo 

o ~ "T ~I 7” 



7CJ* - cj; - 300 



r __ 3Qo f oMU)?-, +'£2 cj n- !°J 

„ ,2 » \3 _ 



7W„ -«K - 300 



'q il 

To determine the stability region, the criterion is 



(7-13) 



Co 3^0 and C, > O u ’ 14) 

Define the dominant roots by a = > f rom Table 2-7 of Chapter II y 

the transformed coefficients of the reduced characteristic equation are: 

P, — C, - ZQ_ 
p 0 — Co-c.a. +a 

The criterion for dominant roots by this definition of dominant roots is: 

P, and j> o o ( 7 - 15 ) 

Computations for ^ =0.3 and ^ = 0.5 from Equations (7-13) and (7-12) 

are shown in Table 7-2 and Table 7-3 respectively. 
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cu n 


C 0 


P 


0.2 


9.13 


0.001217 


0.4 


8.22 


0.0044 


0.6 


7.25 


0.0084 


0.8 


6.34 


0.0135 


1.0 


5.28 


0.0176 


1.2 


4.16 


0.02 


1.4 


2.97 


0.0193 


1.6 


1.71 


0.0146 


1.8 


0.397 


0.0043 


1.9 


-0.312 


-0.0038 



Table 7-2 ^ 



U) n 


Co 


P 


0.6 


5.84 


0.007 


0.8 


4.46 


0.0095 


1.0 


3.07 


0.01022 


1.1 


2.24 


0.00904 


1.2 


1.56 


0.00748 


1.3 


0.87 


0.0049 


1.4 


0.1925 


0.00125 


1.5 


-0.48 


-0.0036 



Table 7-3 f = 



c, 


cX 


3>/ 


J> 

O 


7 


0.0046 


+ 


+ 


6.88 


0.0101 


+ 


+ 


6.81 


0.0166 


+ 


+ 


6.49 


0.0235 


+ 


+ 


6.41 


0.0313 


+ 


+ 


6.27 


0.0397 


+ 


+ 


6.21 


0.049 


+ 


+ 


6.07 


0.057 


+ 


+ 


8.86 


0.097 


+ 


+ 


4.73 
0.3 for 


Example 7-2 


+ 




c, 


<X 


3>, 




6.4 


0.019 


+ 


+ 


6.2 


0.0191 


+ 


+ 


6.01 


0.0216 


+ 


+ 


6.02 


0.0300 


+ 


- 


6.05 


0.0322 


+ 


- 


5.75 


0.036 


+ 


- 


5.61 


0.0375 


+ 


- 


5.46 


0.0383 


+ 


- 



0.5 for Example 7-2 



From the criteria (7-14) and (7-15), the stability and dominant region 



for f = 0.3 and ? = 0.5 are determined from the computation tables. 

Dominant ( O. - f ■ J 
CJ n < 1.8 
tj n c I, OS 



f = 0.3 
^ = 0.5 



Stability 
U; n < 1.8 



Since this computation is a progressing process*, when the stability 
and dominant root region have been located* the compensator parameters 
also have been calculated already. The computation tables not only in- 
clude the compensator parameters, but also the secondary roots. The 
coefficients (D's) of the transformed equation are shown signs only in 
the compensation tables, because in this case only signs are needed to 
determing the dominant region. 
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Assume the dominant roots are chosen as ^ = 0.3, QJ n = 1.0, then 

from Table 7-2, C q = 5.28, = 6.41, p = 0.0176, OK = 0.0313 and 

p 

£ — — — = 0.563. The reduced characteristic equation is 

s 2 + 6.41S + 5.28 = 0 

Then the secondary roots are S = -5.44, -0.97. The closed loop pole- 
zero configuration is shown in Fig. 7-4(a). In Fig. 7-4(a), all the 
secondary roots are to the left of the dominant roots. 

Consider Table 7-2 and 7-3, they indicate 1, this implies the 

compensator is a lag network. For a lag network the zero is close to 
the origin. Then a root to the right of the dominant roots is allowed. 

This dipole near the origin has negligible effect on the transient response 
if the root is close enough to the zero. In this sense, the definition of 
the dominant roots may be adjusted as such that all the secondary roots 
must be to the left of the dominant roots except the root which is close 
to the zero. 

In this example: if f = 0.5, U) n = 1.2 is chosen as the dominant 

roots, from Table 7-3, C q = 1.56, = 7.25 jb = 0.00748, - 0.035 

p 

and Z = •— = 0.214. The reduced characteristic equation is 

s 2 + 7.25S + 1.56 = 0 

Then the secondary roots are -7.03 and -0.225. The closed loop pole- 
zero configuration is shown in Fig. 7-4(b). Compare the pole-zero con- 
figurations of Fig. 7-4(a) and Fig. 7-4(b), it can be easily seen that 
the scheme of Fig. 7-4(b) has better performance than that of scheme 
Fig. 7-4(a). 
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/co 



Fig. 7-4(a) Closed Loop Pole-zero Configuration for the Dominant 

Roots ^ =0.3, aJ n = 1 of 
Example 7-2 



\ 




JCO 



Fig. 7-4(b) Closed Loop Pole-zero Configuration for a Dipole ( 
Near the Origin of Example 7-2. Dominant Roots 
are P =0.5, = 1.2. 



i n o 



7-2 Stability Boundary Lines. 

It has been pointed out that the stability boundary lines for a 
cascade compensator are the same boundary lines of relocation zones in 
Reference 7. In the single section zone if the compensator is to be 
realized by an R-C circuit with pole-zero on negative real axis* then 
the necessary and sufficient condition in all the roots are within that 
zone. That is in the single section zone the compensator, Q =. — — ^ 

c U + PJ 

always can be realized by an R-C circuit for any choice of the roots pro- 
viding those roots are within that zone. In the multiple zone, the roots 
being in that zone is a necessary condition but not a sufficient condi- 
tion. 



For higher order systems there may be several zones for a single 
section compensator. But for dominance, only the zone which is the 
closest to the j Cx) axis is of interest. It has been proved that putting 
a pair of complex roots within the lag area and close to the stability 
boundary lines (uncompensated root-loci) implies a closed loop dipole 
close to the origin. By the same reason, putting a pair of complex roots 
within the lead area and close to the boundary line (pseudo-root loci) 
implies all secondary roots are far away from the dominant roots. The 
question of how far away and how close they are, the dominant root bound- 
ary lines of the definition a = ^ U) n serves a guide for both condi- 
tions. 



7-3 Double Section Cascade Compensator. 

If the dominant roots are expected to be within the multiple sec- 
tion zone, at least two sections are necessary. In general, the design 
is carried out by two steps. First choose an intermediate dominant root 
which is close to the stability boundary lines. Form the transfer function 
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obtained by the intermediate dominant roots, then repeat the process. 

The working formulas for the open loop transfer function with one zero 

in it have been derived in Appendix 2 for an original nth order system. 

The design by intermediate dominant roots has two advantages. First the 

quadratic term as discussed in Chapter 3, equations (3-3) is avoided. 

Secondly the compensator can be designed as a double lead section, a 

double lag section or a lead-lag network. The intermediate dominant 

roots should be chosen close to the final dominant roots because the 

regions which are close to the dividing lines (both uncompensated root- 

loci and the pseudo root-loci) are dominant root regions. 

4-00 

Example 7-3: Given C.J — ■ " m ^~ > the system is expected to operate 

on a pair of dominant roots ^ =0.7, 0O n = 20, and is not to be re- 
duced. 

From Example 7-1 and Fig. 7-2, = 20 is on the stability boundary 

line, then double section is needed. From Table 7-1, choose the inter- 
mediate dominant roots ^= 0.5, U)^ = 22, which is close to the final 

dominant roots ^ =0.7, U>^ = 20. The transfer function becomes 

__ ^8Q(S + !8S) (7-16) 

^ ~ S X C S + 242.) 

From Appendix 2 equation A-2-4, the four equations are 

p- C, = 

8,p + -fa* - Co - if «*.<:, = W,'- 6, ->--§7 

6, p -h 0 o o( - 2 ^ iO t} Co ~ C - 1 — ^ 

6 c p -w,Xo = O 



(7-17-1) 

(7-17-2) 

(7-17-3) 

(7-17-4) 
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From A-2-6: 



z C^- +if U)*] 

r — — (/-lo; 

0 [8.- 8,0)' -fifui,, 1 ] t [-ij’g.f 

where B q , and B£ are the coefficients of the characteristic equation 
of the open loop transfer function (7-16), (-C) is the open loop zero . 

The stability boundary lines are shown in Fig. 7-5a by the rules of 
reference 7. From Fig. 7-5, it indicates that the roots of ^ = 0.7, 

LO n = 20, are within the single lead section zone and close to the pseudo 
root-loci. Then they are expected to be dominant. Substitute C = 185, B q 
97500, Bj = 5280, B 2 = 242, ^ = 0.7, = 20 into equation (7-18), one 

obtains 



From (7-17-4) 
From (7-17-1) 
From (7-17-3) 
Then 



C = 6340 
o 

p = 26 

C L = 240 
<X = 1.4 
Z = 18.6 



The secondary roots are determined from the equation 

14-03 f (oi4-0 — O 



( 7 - 19 ) 



The solution of (7-19) is -210, -30. 

The compensator for the second stage is 




(Sf 

( S "f -2 (e>J 



The overall compensator is 

_ ^ t '1^1 

(s + cs+^&j 

The closed loop pole-zero configuration is shown in Fig. 7-5b. From 
Fig. 7-5b it can be visualized that the response is speeded up by the 
two zeros while the overshoot is small because of the approximately equal 
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jto 




Fig. 7-5(b) Final Closed Loop Pole-zero Configuration of 

Example 7-3 
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distances of the zeros and the dominant roots to the origin. 



7-4 Mixed Compensators. 

For a single section cascade or feedback compensator it has been 
shown that the adjustable parameters are linearly dependent (Matrices 6-9 
and 7-5). In general a compensator has the transfer function 



dependent for either a cascade or a feedback path. The solutions given 
a set of roots are always real numbers. If there are more than one compen- 
sator and in different loops, the adjustable parameters are not linearly 
dependent. The solutions for a given set of roots may be complex numbers 
and consequently the compensators are not physically realizable. However^ 
the subset of roots which gives the solution of real numbers can be found 
if such a subset exists in the stable root region. Note the coefficients 
of the reduced characteristic equations are always real numbers for any 
choice of the roots and for any type of compensators. The reason for this 
fact is due to the complex conjugate roots. Because of this fact the stable 
root region and dominant root region are always defined if they exist. 

The physically adjustable parameters also appear as real numbers in some 
transformations for any choice of the roots, but the individual parameters 
may not be real numbers. Consider Fig. 7-6 which consists of a single 
section compensator in cascade and another in the feedback path, <2f is 
the given plant and , p, , p l and ^ are adjustable parameters. 

From a signal flow graph, the characteristic equation is 




( 7 - 20 ) 



The coefficients, , (X, , -- Q p and b 0 , - - 



b ^ are linearly 
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( 7 - 21 ) 



I 




Fig. 7-6 Mixed Compensators 




Fig.7-6(a) Mixed Compensation of Example 7-4 



I 



134 



Consider the first term of (7-21), it can be expressed as follows: 

(i+P«)(S + P*)Q = [S*+ C P.+ PJS t (7-22) 

Let 

P 1 + P 2 = b (7-23) 

P, = 0_ (7-24) 

where b and are always real numbers for any choice of roots of (7-21). 
However, p and p are not necessarily real. To find the real solution 
of p ; and , the roots must be subject to some restrictions. In this 

case it is 

b > (7-25) 

The inequality (7-25) defines the realizable region. For a passive net- 
work, additional constraints are necessary, that is 

b a >- a (7_26) 

The same constraints of the roots can be applied to other cases. 

Example 7-4. Given (j — — , The velocity error 

H £(S+JrJ(S + lo) z 

constant is not to be reduced. The band width is less than £o rad/sec. 
From the specifications dominant roots of € =0.5 Cj = 10 are chosen. 

For cascade compensation alone, more than two lead sections are required 
by the Ross Warren method 7 . A scheme of Fig. 7-6 (a) is attempted. Be- 
cause of the gain requirement, only second derivative feedback is used 
and the parameters o( , £ and p are related by the following relation 

C{ = ~ (7-2 7) 

Z 

There are three free parameters ( , p , ), three roots are art it - 

rary. The constant controlled characteristic equation is 

s £ '+ (2<r+pJ5 4 + (aoo + ^oop + \o*L.)£ V (£ooi2vof> tto%f>)s 2 
_|_ (iroop + \<f*) 5 -f- I0 S p — O (7-28) 
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Let (-S^), (-S 2 ) and (-S^) be the arbitrary roots and define the 
following terms for simplification. 

F, = 5,+5^S J C 7 - 29 ’ 1 ) 

F 2 = _S,S t + S 2 Sj + SjS t ( 7 - 29 - 2 ) 

F 3 = $,S^ ( 7 - 29 - 3 ) 

The root-coefficient relations are: 

i TOO 4 p =. F, 4 C, (7-30-1) 

ZOO 4 i'OOp -t I0 S '^= C* + C,F, + Fi (7-30-2) 

25 + ^oop 4(o r &j3 = CoF, +C,F i + F, (7-30-3) 

Z5 p 4*10^ = Co F 2 + C, F, (7-30-4) 

10 * P = C„F 5 (7-30-5) 

Where C, and C are the coefficients of the reduced characteristic equation 
1 o 

which defines the constrained roots (-S^) and (-S^) by the equation 

i> Z + C,S f Co — o ( 7 - 31 ) 

From equations (7-30-1), (7-30-2), (7-30-3) and (7-30-5), C is 

o 

obtained : 

Fj ( F.f, -2SFj *10^ cl - lO r [F, , Fj -irF.Fj + '°‘> ;]Co 

i I0'° [ f, f A - 2$ F z - Fj -t soo] = O (7 ' 32) 

Assume (-S^) and (-S 2 ) are the dominant roots and denoted by ^ and CJ n # 
(-S^) is arbitrary and must be greater than • Assume S 3 = ZfU n 

For ^ =0.5, LaJ^ = 10 from (7-29) 

F, = 20 

h = 200 
= 1000 
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Substitute into (7-32) one obtains, 



C = + 2080, - 31.6 

o 



Take 
From (7-30-5) 
From (7-30-1) 
From (7-30-2) 
From (7-30-4) 



C = 2080 
o 



P = 20.8 



C x = 25.8 



10^= 2076 
OK, = 4.32 

The reduced characteristic equation is 

5 2 + ^S8S + zo8o = O 

The constrained roots are -S. , -S, = -12.9 + j 43.5 

4 5 — 



The five roots are: 



■S , -S, = S + j 8.66 



-S 



= -10 



The compensator is 



3 

-S 4 , -S 5 = -12.9 + j 43.5 



_ 48Z ( S-K8iJ 

( - T C~~ S + 2.0,8 



L — o.odcna 



The closed loop pole-zero configuration is shown in Fig, 7-7. For other 
choices of the arbitrary roots, better results can be obtained. 




Fig. 7-7 Closed Loop Pole-zero Configuration 
of Example 7-4. 
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VIII. CONCLUSIONS 



8-1 Summary of Results. 

A new method for design of control systems in the time domain has 
been described. The method is applicable to any system configuration and 
is able to handle any number of variable parameters. The design procedure 
is based upon a set of simultaneous algebraic equations, and is well suited 
to implementation by digital computer. The method is flexible and has more 
freedom in design than other methods. The values of adjustable parameters 
are computed from analytical expressions, so any desired degree of accuracy 
can be obtained. 

The main advantage of this method over those methods available in the 
literature lies in the fact that the method is a direct syntheses from the 
closed loop pole-zero configuration. No trial and error is needed. The 
specifications and restrictions of the compensator confine the arbitrary 
roots within a region on the S-plane, once this region or regions have been 
determined, the rest of the design work is just a computation. The compen- 
sator obtained by this computation guarantee that the specifications and re- 
strictions are satisfied. 

The main difference between this method and other available methods 
is that the domain of the system can be any number of variables. Theoreti- 
cally, the more variables available, the more design freedom one has. In 
the root-locus method, only one parameter is varied others must be fixed. 

In Mitrovic's method, only two parameters are varied. This method is able 
to handle any number of variable parameters simultaneously. 
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8-2 General Description of the Method. 

The design procedure can be summarized as shown in Fig. 8-1. The 
analytical expressions of the dependent variables of commonly used com- 
pensators have been derived in the previous chapters. In general for the 
case in which coefficients of the characteristic equations as variables^ 
the dependent variables are always in linear relations one to another. 

A compensator of the transfer function as (8-1) is in the forward path or 
in the feedback path alone, the coefficients a’s and b’s are always linear 



relations in the partitioned formulas. Those linear algebraic equations 
are always consistant and therefore always have a solution. For the case 
in which compensators are in both forward and feedback paths, a quadratic 
expression is always obtained. But this quadratic expression does not 
imply the dominant and realizable regions are smaller. It shows only the 
natural algebraic relations. For example, the coefficients of the reduc- 
ed characteristic are always real no matter how complicated their functions 
are. 

8-3 Extensions. 

The method can be extended to sampled data systems and adaptive sys- 
tems. The modifications are the different interpretations of the variables. 
For example, in the case of sampled data systems, if £ transform is used 
then the criteria of stability, dominance and realizability are different. 
However, the general idea of approach can be applied. 

In the previous chapters, the arbitrary variables are chosen as the 
roots of the system. This choice is not necessary. Theoretically, any 




(8-1) 
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Pig. 8-1 General Description* of the Method 
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combination of "r n variables can be chosen as the arbitrary variables such 
as a combination of roots and variable parameters. 

This method can be applied to any mode of linear operation in dis- 
continuous operation. By proper choice of the compensator, the dominant 
roots of two modes can be chosen in a relation as the designer wants. 
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APPENDIX I 



CONSISTENCY 6c LINEAR DEPENDENCE OF DERIVATIVE FEEDBACK 
1. Introduction, 

The coefficients of the characteristic equation of a linear control 
system are functions of the system parameters. If all the coefficients 
can be adjusted by a set of the parameters, then the roots of the system 
can be controlled. Conversely, if the roots are chosen, then the system 
parameters are forced to have certain values. The conditions of using 
system parameters to satisfy the coefficients of a characteristic equation 
for a given set of roots depend upon the number of adjustable parameters 
and the consistency of the simultaneous equations. Theoretically, if all 
the derivatives of a system can be obtained from the output of the system, 
then by feedback as shown in Fig. A-l-1, the system roots can be chosen 
arbitrarily providing no restrictions are placed on the forward gain K 
and the control parameters h's. This can be shown as follows. The charact- 
eristic equation of the system with feedback is: 



Assume the roots of the system are chosen from the considerations of 
the specifications, and the required characteristic equation is formulated 
as equation (A-l-2) 



For the system to have the chosen roots, then the corresponding coeffi- 
cients of equations (A-l-1) and (A-l-2) must be equal. Then 






2 . 



+ f\i 5 -h A,S> i f\ 0 ~ O (A-l-2) 
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Fig. A-l-1 Block Diagram by Feedback of Derivatives . Gr is 
Transfer Function of the Plant ^ K and 
are the System Gains. 



) 
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, . < 



k,K a - 1\ 0 



kzh, + ^ I (A-l-3) 

Ki ^ o-i ^ ^-n -1 A n- • 

Equation (A-l-3) always has a unique solution if there are no re- 
strictions on K f s and h's. And in the particular case, each parameter 
h controls one coefficient independently and the forward gain controls 
the constant co-efficient. If the control parameters have certain restric- 
tions, then the system roots cannot be arbitrarily chosen, but have to be 
subject to these restrictions. For example, if is specified from 

the consideration of steady state accuracy, then the product of the roots 

must equal to K,K 0 . If h cannot be obtained (h =0), then the sum 
1 Z n-1 n-1 

of roots must be equal to a^ . The restrictions on the control para- 
meters restrict the roots in a certain constraint. The more restrictions, 
the higher is the degree of the constraints of the roots. In a practical 
system, the higher derivatives are very noisy and very difficult to obtain. 
This is a common restriction for derivative feedback. 

In Fig. A-l-1, it was assumed that there were no zeros in the forward 
path. If one zero is in the forward path, then the allowed highest order 
of the feedback derivative is (n-2). In general, if >n" zeros are in the 
forward path, only the derivatives of the order up to n-(M+l) are allowed 
to feedback. This can be proved as follows 

Refer to Fig. A-l-1, assume m zeros in (jp , and the highest deriva- 
tive has an order of p . Then the characteristic equation is 

5 n + (a„., +h rtln ) s' 1 "' + - - - =<9 (a-1-4) 
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The subscript ( f\\rs ) of h r+m indicates the highest order of the feed- 
back derivative and the number of zeros. From equation (A-l-4) 

r + hi c n — i 
r < n - ( m + ij 

Therefore if ni zeros are in the forward path, then only n-(m+f ) 
control parameters can be introduced, if the feedback path is as shown in 
Fig. A-l-1. However, other schemes of lower order derivative feedback are 
possible to introduce (n-1) control parameters and consequently control all 
the coefficients of the characteristic equation. The following sections 
are discussions about these possibilities. In the following sections, two 
terminologies are defined as follows: 

Control variable is the variable of the system which can be measured. 

Control parameter is the adjustable parameter in the system. 
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2. Feedback of the First and Second Derivative from the Output to 
Several Points in the Forward Path. 

Since the first and the second derivative are usually allowed to 
feedback for high order systems and can be obtained, only these two deriva- 
tives are considered. The object of this discussion is to adjust all the 
coefficients of a nth order characteristic equation to desired values. In 
order to accomplish this, n control parameters are needed. As the forward 
gain K is considered as one control parameter, (n-1) control parameters 
from the feedback path are necessary for this particular scheme. Using 
first and second derivative feedback around the plant intorduces only two 
control parameters. However, if there exist junction points between the 
energy-storage elements in the forward path, such that signals can be fed 
in or picked up, then additional control parameters can be introduced by 
feeding derivatives to those points. As two control parameters have been 
introduced by feeding derivatives to the error channel, then the neces- 
sary number of points between the energy storage elements in the for- 
n 3 n 2 

ward path is — ^ — (if n is cdd and — if n is even) since two 
derivatives can be fed into each node. Briefly, this can be summarized in 
the following two statements. 

n- 3 

(1)* If a nth order system has at least — j- available points 
n-2 

(if n is odd, and if n is even) between the energy storage elements 

in the forward path, then all the coefficients of the characteristic equa- 
tion can be adjusted to arbitrary values by introducing first and second 
derivatives from the output to each point providing the two derivatives 
are defined and the feedback paths are proper. 

★all statements in this discussion may be considered as rules. All of them 
are necessary conditions, but not sufficient conditions. 
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(2) If a nth order system has at least n-2 available points 
between the energy-storage elements in the forward path, then all the co- 
efficients of the characteristic equation can be adjusted to arbitrary 
values by introducing only first derivative feedback from the output to 
each point, providing the feedback paths are proper. 



^ I (,) -I (2> 

L k h- ^ l k 

{ M 

I zn. product of the kth possible combination of non- 

K 



3. Formulation of Characteristic Equation for Pure Derivative Feed- 
back, 

From the signal flow graph, the determinant of a system can be 
obtained by the following rule, 

A = I 

( 'U 

where 

touching loops 

I j = loop transmission 

For this particular derivative feedback system^ two terminologies 
are defined as follows: 

Loop gain = numerator of each loop transmission. 

In Fig, A-l-2, Loop gain of L^ = h^P 3 

Loop gain of L 2 = h 2-2 P 3 
Loop gain of L 3 = h 3 P 1 P 2 P 3 
Loop gain of L^ =* P 1 P 2 P 3 

Characteristic function = denominator of each of the forward transfer 
functions, or the product of them. 

Consider Fig. A-l-2; from the graph rule 



.A 



| + L. + L 2 + L > + L. 



= I + 



P 3 b, 



Q. 



4- i 



4 PPz P-3^ J. P, Pz P, 



u ' Q'QA 



Q&Qj 

Q,Q 2 Q 3 + Q,Q* ( Pjh.J 4 Q.CP.P,^) 

+ 4 Pi PiP3 



(A-l- 5 ) 



The first term Q ( 4 last term = characteristic equation 

of the original system. 
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Fig. A-l-2 Formulation of the Determinant 

Q,(s) h;(s) ore poly no m « na I fp S 



The second term = loop gain X nontouching chara. function. 

The third term = loop gain X nontouching chara. function. 

The fourth term = loop gain X nontouching chara, function. 

From equation (A-l-5), the formulation of the characteristic equa- 
tion for pure derivative feedback may be stated in the following rule. 
The characteristic equation of a pure derivative feedback system - 



(0 

K 



(z) 
^ K 



loop gain X 
loop gain X 



characteristic equation of the original system + 
characteristic function of nontouching loop + 
characteristic function on nontouching loop. 

Example A-l-1: Devise a derivative feedback scheme of Fig, A-l-3, 

such that all the coefficients of the characteristic equation can be ad- 
justed. Assume the required characteristic equation from the specifica- 
tions is 

/i 4 s* + a 3 fl.s + ft 0 =- o 

find the relations of the control parameters and the coefficients of the 
given characteristic equation. 

Since intermediate points are available, then n-2 = 5-2 = 3. From 
rule (2), it is not possible to use first derivative only. 

By rule (1) = 1. Then a combination of first and second 

derivative can meet the requirement. There are many possible ways to 
accomplish this since only 4 control parameters are needed. 

Two schemes are listed as follows: 



(a) 


H i = V 


^2 ” ^2 S + ^3 S 


H 3 - V 


(b) 


H i - V 


H 2 ' h 2 s2 


h 3 = h 3 s2 + V 



Assume scheme (a) is chosen. The chara. equation is formulated in 
a tabulated form for convenience. 
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1 




// 






X 




t / 




Fig. A-l-3. Example A-l-1 - Mark "X" denotes signals are allowed 

to feed in or pick up 
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Required Chara, Equation 



Original Chara. Equation 
Gain x Chara. Nontouching 

\S 2 K, K l Kj (s+<*) 

Gain L3 x Chara. Nontouching 

(h i s l +h J S)K 2 ^*(S+Pi) 



Gain x Chara. Func. Nontot: 

h,s Kj *(s+r,)(s+r 2 )(s+frj 



r 4 


3 

| 


! 


H 




1 

zrj; 


i 

^Pc^-Pk 


i 

rr P; 1 

/V 








-»W$ 




k.Wu 


tXK 3 °< h + 






k 2 M. 


K*K}P| 

+ l<2^hj 




; KjH* , 


KjMpx 


MP,P« + 


Wi*). 



* WX 



+ iy f \w s r,)h, 



The equations for the coefficients are 

K)l7, 

2- Pc Pj + K, Ki t'i +■ Kj h 2 + ^3 f fi + Pi F Pj} h , = ^3 

^ U.C. ^ K.KVKs^h^. + 1C* Kj P K; hj + k} (P,f t + Pifs i fl. 

TF P- <- K 0 K, K, K s +• k,r, P, -/!, + Kj PJJj h, — A \ 

ot. K 0 fC, K 2 fCj = 



The first four equations of (A-l-6) can be written in matrix form: 

-ZP- 



fc 3 000 



KjU'lWi) K,fc, O IC.K.Kj 

K z (C 3 p, K z K 3 

Ks Pi K Fj, Q U L - p Q 



r “7 

fl 



*lz 

Hi 

»u 



Al'ZrM 



Pi P; Pi 

\5L 



In (A-l-7) whenever the determinant is not zero, there always exists 



a solution of the control parameters h’s for arbitrary values of A. 

In the above example, the table of the characteristic equations gives 
a clear picture of how the coefficients are controlled by h-parameters. 

Also equation (A-l-7) shows that the h-parameters are linearly dependent 
to the coefficients of the given characteristic equation., It is this 
linear dependence which makes the arbitrary adjustment possible. In 
equation (A-l-6), the constant term is controlled by the forward gain only. 
This is the nature of derivative feedback, because the derivatives control 
the coefficient of the terms of a differential equation. Any derivative 
feedback from the output has the scheme of properties of linear dependence 
of h-parameters and the independence of forward gain K. 

Example (A-l-2). Fig. A-l-4 is a third order system, one point between 
the energy-storage elements in the forward path is allowed to feed in signal, 
then the system roots may be arbitrarily chosen. Assume the roots to be 
-3 and -0.5 + j 1. K^, h^, h^ can be determined. The characteristic equa- 

tion for the given roots is; 



S 



3 



+ 4S 2 + 



4.25S + 3.75 = 0 



Main loop 



1 



1 



0 



K 



1 



L^: loop gain = 



chara. function 
not enclosed = 1 




1 * 2 : loop gain = h^s 
Chara. function 
not enclosed = s + 1 



h 



1 



h 



1 



153 




Fig. A-I-4 Examjtrle A-I-2 



X 



tf 



"O — O" 



ii — 





SO 


i 


J 




^ — V 



o 



I — 



■ 2Q - 

sc-s+a; 




Fig. A-I-5 Example A-I-3 
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Set the corresponding coefficients equal, variable parameters are solved. 

= 3.75, h 1 =3, h 2 = 0.334 

Example A- 1-3. In Fig. A- 1-5, it is assumed that signal can be fed in at 
the point "X" in the forward path, and the forward gain of the system is 
1500 based on the steady state specification. According to rule 1, = 1, 

then by first and second derivative feedback, the system roots can be chosen 
arbitrarily providing the product of the roots is equal to 1500. 

The uncompensated system is unstable. The feedback paths are chosen as 
shown, and the roots are assumed to be S = (-3 + j5), -6, and -7.36. 

The characteristic equation for the chosen roots is 

S 4 + 19.36S 3 + 158. 5 2 + 720S + 1500 = 0 

The system coefficients 



Uncompensated system = 1 


10 


31 


30 


1500 


L 1 : ( hfS+^S 1 ) 1500 




1500h 2 


1500h 1 




L 2 : hjS 30(5 + *Xi+$-J 


30h 3 


240h 3 


450h 3 





The coefficients relations are: 

30h 3 + 10 = 19.36 

1500h 2 + 240h 3 + 31 = 158 

450h 3 + ISOOhj + 30 = 720 

Solve equation (A-l-8), get (A-l-8) 

h = 0.367 
h 2 = 0.0347 
h 3 = 0.312 
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5. Feedback from Intermediate Points. 

It can be shown that if the output of each block in the diagram (A- 1-6) 
is chosen as the state variable of the system, then there exists a linear 
transforamtion between the state space and the phase space (output and all 
its derivatives) of a system. In Fig. A-l-6, the state variables are 
chosen as shown, the relationship between state space and phase space is 




Fig. A-l-6 State Variables 



derived as follows: 

© c 



(A-l-9-1) 



X, J_ v __ C _ (A-l-9-2) 

TT = 3; “ ** 



Kx Ko 
“i+b J 



+ b 



*■> ~ kJ 



(<£) 


, -L f49 \ 


( cttj 


+ Ku l city 



(A-l-9-3) 



Equations (A-l-9) are the required relations. It shows that the out- 
put of each block is a linear combination of the output Q and its de- 
rivatives. 



156 



With this concept, therefore, the feedback from intermediate points 
is equivalent to derivative feedback from the output. Fig, A-l-7(a) is 
equivalent to Fig. A-l-7(b). 

This equivalence also can be shown by block diagram manipulations. 
Because of this equivalence, if the signals flowing in the intermediate 
points can be measured, then those signals can be taken as control vari- 
ables and the analysis is the same as for derivative feedback from the 
output. Moreover, if the derivatives of those signals are defined and 
can be obtained, they also can be taken as control variables. 

However, if consistency and linear dependence are considered, the 
feedback paths and the number of control parameters have many restrictions. 

(1) Consider the scheme of Fig. A-l-8, it introduces sufficient 
control parameters but the coefficients of the characteristic equation are 
not linear combinations of the h-parameters . It consists of two indepen- 
dent nontouching loops, the contribution of these two loops to the whole 
system is the product of the two loop transmissions. So the coefficients 
are quadratic functions of the h- parameters. This can be shown by formu- 
lating the characteristic equation as follows: 



s 3 


s 2 


s 1 


s° 


Original system 1 


cub 


ab 


K*K b 


Lj ; K b h,S (S+aJ 


K b h, 


KPh, 




^2 Kr (S 4- bS) 


KA 






L 1 X L 2 




WA 





If the required characteristic equation is 



S 3 + A 0 S 2 + A,S + A = 0 
2 1 o 
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Fig. A-l-7(a) Feedback from Intermediate points. 



tCo. 

5(S+cO(S + bJ 





Fig. A-l-7(b) Equivalenence of Fig. A-l-7(a) 

) 
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\\ 

X 




Fig. A-l-8 Quardrtic Function of h parameters 



\\ 





Fig. A-l-9 Example A-l-4 
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then 



Kuh, 4 Kab t = fl 2 - <0+ bj 



K a 0bi 4 ^,bfl 2 + (A-l-10) 

K a k b = /!„ 

The solution of equation (A-l-10) for real values of and ^ implies 
the restriction of the roots in a certain region. In general, if there 
are " /[ " independent feedback loops in a system, the coefficients of 
the characteristic equation consist of the terms h^ 9 (h^h^,) s . . . (h^. . .h^) . 

This leads to the following statement. 

If all the coefficients of the characteristic equation of a derivative 
controlled system are a linear function of the control parameter, then all 
feedback paths must be interconnected, (or touching). The above state- 
ment may be stated in another way: If all the coefficients of the character- 

istic equation of a derivative controlled system be a linear function of 
the control parameters, then those points from which the signals are 
taken as control variables cannot be used as feed-in nodes. 

Example A-l-4. In Fig. A-l-9, assume the signal flowing at the point 
marked "X" can be measured, then it can be picked up as a control variable. 
Assume the feedback paths are as shown, they are interconnected. The 
characteristic equations are: 

^ ? 

A s 0 
o 

K «*b 



3 

Required characteristic equation = S 


4 A 2 S 2 


+ 


Characteristic equation j 


a-t-b 


ab 


L 1 - h,S K a K„ 




U„h, 


L 2 “ U 2 + bs) 


K«h> 


K ahx b 



The coefficient relations are: 

Ko = A z - (Q+ bj 

K a K b h, + k Q bh z = 

Ka Kj, — Ao 

and can be solved for arbitrary value of A's, 



(A-l-11) 
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(2) If zeros are involved in the forward path, the feedback 
from the intermediate nodes may introduce the problem of consistency „ 
Consider Fig, A-l-10, the signal at "X" can be measured and then can be 
taken as control variable. The characteristic equations are: 



Required Chara. Equation 


s 3 


2 

+ a 2 s^ 


+ A^ S + A = 0 

1 o 


Original loop 


l 


a + b 


ab+k k. K K. 

a b a b 


l 2 * b l SK C{ K b (S+^) 




K<i Kb hj 


hz 


Lj 5 CS-hoiJ(S Z fbsJh, 


l^h, 


fc a (s+bjh ( 


Kqdb h, 



The coefficient relations are: 

M, = i 

K« h, = /!*- i4 + bj 12 

K a oCb h, + “ A t 

KaK h oC = Ao 

Equation (A-l-12) are linear, but are not consistent for arbitrary 
A's, since the number of equations is one more than the number of unknowns 
(h^ and t^). Substitute h^ from the first equation into the second and 
third equation, obtain: 

+ (^ + bj = A 2 -C^ + bj 
K a K b °L K + <* b = a, - (ab + Kq K b ) 

Manipulate, get 

ckA z - A, — oi (o{ + Q - 3bJ - ab -■ (A-l-13) 

Equation (A-l-13) is the restriction of the roots for this systm, i,e. s if 
equation (A-l-12) has a solution, then equation (A-l-13) must be satisfied. 
Therefore, when a system has the required number of control parameters and 
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"X" 




Fig. A-l-10. Effect of Zeros to Feedback 




Fig. A- 1-11. Feedback of the output to more than one point. 
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the linear relationship of the coefficients, the choice of roots may still 
be subject to some restriction if the feedback is not proper. In general it 
can be shown that if the loop transmission has the same order (s) for both 
numerator and denominator, then it introduces the restriction of the roots. 
This leads to the following statement. 

If the roots of derivative controlled systems are to be chosen arbi- 
trarily, then for each loop transmission, the order of the numerator must 
be at least one less than the denominator. 

(3) There is another obvious inconsistency which arises from 
feedback of the output to more than one point. Therefore, if an inter- 
mediate signal which is not followed by pure integrator, be fed to an- 
other point, the equations of the coefficients is inconsistent for arbi- 
trary choice of roots. Fig. A-l-ll is an example. 

The characteristic equation of Fig. A-l-ll is 



original loop 

L x : h ,K,(S+oJ(S+b) 

V hjk.KgS 



S' 

/ 



Qib+c ab+bC+ba K,K Z & +<*bC 



M, 



{c\ibjK,h> 



abk, h i 



K,K z hi 

There are three equations but two unknowns. This leads to a con- 
straint of the roots. 

(4) Since the feedback from an intermediate node is equivalent 
to derivative feedback from the output, then it is possible that the para- 
meters introduced by both control variables may have the same effect on 
the system; i.e., the parameters may not be independent. This can be 
easily detected by just looking at the loops. In Fig. A-l-12, h^ 
has the same effect as h^ and h^, since h^ is a combination of the first 
and second derivative of the output. That is, h^ enters the coefficient 
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< 




of the first and second derivative terms as h^ and h^. Therefore in 
Fig. A-l-12, the system actually has only two control parameters. 

However, if the derivative of the signal at "X" can be obtained 
and the scheme changes to Fig. A-l-13, then each parameter plays its own 
role. 
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) 




Fig. A-l-12. Two parameters have the same effect. 




Fig. A-l-13. Each parameter has it's own effect 



) 
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5. Summary. 

Derivative control is an effective method to control the coefficients 
of the characteristic equation of a system, since it uses derivatives of 
the system variables to control the coefficients of the characteristic 
equation. Since an intermediate signal feedback has the same effect as 
derivatives from the output, they are considered as in the category of 
derivatives. The criteria in Section I also applies to the intermediate 
signal feedback because of the condition (1) in Section 4. The inter- 
mediate signals do not introduce more control parameters but increase the 
possible ways of obtaining derivatives and the feedback path. 

In general, one constraint of the roots is almost always required in 
a system, that is, the product of the roots equal to the system forward 
gain, because the gain is usually specified from the stiffness of the 
system. This is only one degree of constraint, it does not introduce 
much difficulty. But if additional constraints are involved such as 
quadratic functions and inconsistent equations, then the solution is not 
so obvious. Therefore, it is logical to look for a set of linear and 
consistent control parameters. From the discussion in the previous sec- 
tions, it may be summarized as follows: 

(1) The number of h parameters must be equal to (n-1) 

(2) h - parameters must appear only in the coefficients from the 
first derivatives to the (n-1) derivative, 

(3) h-parameters must be linear combinations of the given co- 
efficients. 

Example A-l-5. Fig. 2-1-14 is a 6th order system. Assume two nodes 
between the energy-storage element are allowed to feed in or pickup signals, 
and the derivative of the signal at "X'\ can be obtained. From the 
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6-2 

criteria of intermediate nodes, ( — ^ = 2) the system coefficients can 

be adjusted to arbitrary value by first and second derivative feedback. 
Fig, A-l-14 (a) (b) and (c) are the possible proper feedback paths. Fig, 
A-1-I4(c) again has 5 possible combinations of H # s. They are 

H !=h,5+hi \ « h,5( or h,) H x =h ( >+h 2 



H 2 = h 3 s+h f 



H, = h z s + h 



r 3 



H 2 - h 3 5 



H = hj-S ( 



Hj - h 4 S 4 hj-S 



h 3 = h 4 S 



Therefore, there are 8 possible ways to control the coefficients of the 
characteristic equation of this system. The choice of the scheme depends 
upon the specifications and perhaps the engineering judgments. 
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* 

Fig. A-l* 14. Example of A-l-5. 



! 
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APPENDIX II 



PARTITION OF A nTH ORDER SYSTEM WITH ONE ZERO FOR SINGLE SECTION CASCADE 

COMPENSATOR 

Take a 4th order system as shown in Fig. A-2-1. 




Fig. A-2-1. Cascade Compensator 
The controlled characteristic equation is 

S^+ (ft, + pJ5* + (p.tfpj 4 fV S 

+ K‘=<S i + K°<(H c.js +K pc =0 

s i_ + (( 3 J tp; 5 ‘"+(d+^pji 3 +(p,+ pa + ^5" (a . 2 . 1} 

+ (ftp + Kp -i K^C )S + K pC =0 

The characteristic equation of the uncompensated system is 

S^+ p x s'+ ( (3 f KJ-S + KC = o < A-2-2) 

In order to express the coefficients of (A-2-1) in terms of the coefficients 
of the uncompensated system in a simpler form a defines 
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( 



(A-2-3) 



ft = 

ft = B, 

P, + K = B, 

KC = B 0 

Substitute (A-2-3) into (A-2-1) and manipulates, one obtains 

S s + ( B 3 + PJS + + ( B, + pB 3 ) S 3 + ( B, - I 2 + |2=ct 5 s 

(A-2-4) 

+ ( PB, + *B 0 )S + B„p = O 



Set the coefficients of (A-2-4) equal to the corresponding root-coeffici- 
ents and denote the arbitrary roots by ? and U) n 9 one obtains 



(A-2-5) 



Where C q> and are the coefficients of the reduced characteristic equa- 
tion* 

By Cramer* s Rule. 
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_ Bp [s, - 4 - ( I- ^ + Si) ^ 1 ~ T° 

[B.- B^n +*f [-2fB 0 + B,H- 63^+2^ J 



(A-2-6) 



% -« 

<p. =1 

X - 2 f 

% - l-4-f 

170 



Define : 






f°r fC > £ 



For an uncompensated nth order system which has one open loop zero at 

(-C) and has the characteristic equation with coefficients B B 

n-l* n-2’ 

. ...B , then 
o 




g 

[ M ?. * t % + % B 4 w n V] + ig BA) b 4 w> ] 



(A-2-7) 
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